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Current Trends in the American Physical Society 


Recently, there has been considerable agitation in the American Physical Society for 
greater participation in the affairs of the Society by physicists in industry. At the sugges- 
tion of the Council of the Society there is given here a brief review of the discussion of the 
past few months. There follows first a note by Dr. Karl K. Darrow, Secretary of the Society, 
giving the background of this discussion and explaining the Council’s attitude regarding 
the formation of divisions within the Society. Then follow excerpts from letters written 
by members of the Society which indicate the general trend of the discussion. Since the 
letters were not sent in for publication no names are attached to the excerpts. If space 
permits additional comments will be printed as they are received.—The Editor. 


In recent years there has been a great expan- 
sion of activity in the science of physics; the 
number of those who cultivate it has grown con- 
siderably, and the diversity of their specialties 
has increased correspondingly. Moreover a stead- 
ily rising proportion of physicists has been em- 
ployed by industry, a lesser proportion (though 
a rising absolute number) by universities. These 
trends seem likely to go on for years to come, 
perhaps indeed at an ascending rate. It is reason- 
able to inquire whether the structure and the 
procedures of the American Physical Society, 
which have altered but little in the last thirty 
years, may now be modified with advantage. 

This question has been considered by the 
Council of the Society from time to time for 
several years. In 1931 the following article was 
inserted into the Constitution of the Society to 
provide a mechanism for the formation of 
Divisions within the Society: 


“Article IX. 1. The Council may, upon petition by members 
of the Society, form a Division within the Society charged 
with the advancement and diffusion of the knowledge of a 
specified subject or subjects in physics. 2. Each Division 
shall elect an Executive Committee, the Chairman of which 
shall report its activities and needs to the Council. 3. Any 
Division may be dissolved at the discretion of the Council.” 


At each of the Council’s last few meetings, 
including a special meeting held on June fifth, 
there was discussion of the opportunity afforded 
by this Article and the ways of profiting by it. 
These recent discussions were stimulated by 
several letters from members inquiring whether 
the Society could give a greater place and role 
in its activities to physicists in industry and in 
various branches of applied physics. Fortunately, 
in the spring of 1943, advantage was finally taken 
of Article IX by a group organized by Dr. L. 
Marton of Stanford University, who wrote and 
circulated a petition for the establishment of a 
division of Electron and Ion Optics. 
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Signed by forty of the Fellows and Members 
of the Society, this petition was accepted by the 
Council on May first, and a committee was 
appointed to prepare a scheme for getting the 
new Division under way. Further details need 
not be given in this place, for the entire Society 
has already received a description (published in 
the August Journal of Applied Physics) of the 
plans for the new division and of the procedure 
for enrolling therein. This first division will give 
to the Council and the Society a welcome chance 
of observing how a division works and deducing 
how future divisions should be administered. 

The Council will welcome suggestions and 
propositions for other divisions, each to be 
limited in scope to a particular field or fields of 
physics. The limitation should be reasonably 
clear-cut; but no one need apprehend that by 
defining the scope of one division he is cutting 
himself or his colleagues off from other fields, 
for it is intended that any member of the Society 
may enroll in any or all divisions. It is realized 
that in the present conditions, very many 





physicists are entirely absorbed in urgent work 
and cannot give any time to the activities of a 
division. The proverb ‘‘In time of war, prepare 
for peace”’ is nonetheless valid here. We shall be 
designing an adequate framework, in which the 
work of the Society can develop as soon as 
happier times arrive. 

Letters have been received conveying the 
idea of a ‘‘Division of Industrial Physics.” The 
Council is disinclined to favor this idea, which 
incidentally appears to be precluded by the 
language of Article IX. One of the things most 
greatly to be desired is a unification of the 
physicists called ‘‘academic’’ and the physicists 
called “‘industrial.’’ It is important to avert the 
danger of a lack or loss of interest by either 
group in the problems and the enterprises of the 
other. Now, the establishment of a ‘‘Division of 
Industrial Physics’? would work in exactly the 
opposite direction. That way lies the peril of 
forming two contrasting groups, whose aims and 
ambitions ought to be alike but would infallibly 
deviate more and more as the years go on. 

—KarL K. DARROW 





Excerpts from Letters of Members of the American Physical Society: 


Need for Greater 
Unification in Program 


I have heard it said for a 
long time that there is dis- 
satisfaction on the part of 

industrial physicists with the conduct of the affairs of the 
society, but I have yet to meet such individuals personally. 
I don’t mean to imply that they don’t exist, only that I 
have no personal experience with them. 

Personally, I think the dissatisfaction stems largely from 
the fact that industrial physics covers such a wide range of 
topics that it is hard even among the industrial physicists 
to get together much of an audience for any highly spe- 

-cialized paper analogous to the usual contributed papers 
of our program. That is not the fault of the society, but 
simply a statement of the basic fact and I do not see much 
that the society can do to correct it. I think it should 
always be emphasized that the society welcomes papers 
from industrial physicists and welcomes their taking the 

. initiative in setting up special symposia on topics of 
especial interest to them. 

I feel very strongly opposed to any movement to have 
a division of industrial physics or any other such plan 
which would act to emphasize any division of interests of 
those physicists who are employed in industry and those 
employed in academic institutions. We have too much 
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separation now, so what we need is means to bring them 
together instead of means to take them apart. 

It seems to me that the society has done well in fostering 
so successfully the Journal of Applied Physics. It seems 
that the people who are thinking of a separatist movement 
should be encouraged to take the initiative in planning 
program material for the American Physical Society. 





Reciprocal Benefit 
Derived from Contact; 
Sociability Needed 


Personally, I have never 
felt that the Society discrimi- 
nates against industrial phys- 
ics and those who practice it. 
. . . The chief aim of the American Physical Society has 
been, and I believe should continue to be, that of promoting 
the development of physics by bringing physicists together. 
Every physicist knows that the various branches of physics 
are parts of the science as a whole and fundamentally are 
mutually dependent. At this time, more than ever before, 
the Society should use its influence to discourage the prac- 
tice of the various branches of physics by groups of 
“isolationists.”’ 
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It is obvious that I do not believe there should be a 
“Division of Industrial Physics.”’ All physicists learn their 
fundamentals from the academic physicist. The two groups 
should have much in common. The academic physicist 
should be much interested in the progress, success, and 
welfare of the industrial physicist because the latter is his 
chief product. The best way to learn how, and what, to 
teach future physicists is to keep in close touch with those 
out in the industrial field. If the academic physicist is not 
directly interested in carrying his research to practical 
ends he should, and no doubt does, desire to produce men 
who are interested in applying the results of academic 
research to the development of things with which people 
can make life and thé world more livable. 

I suggest that the Society should make the meetings 
more like social gatherings and less like dry old class room 
meetings in which the “professor” is away from his home 
apparatus and with his “‘lesson’’ poorly organized for 
teaching purposes. It might be a good plan for the Society 
to require the presence of every speaker at several 
luncheons and dinners. Instead of one or two tables there 
should be one for each group of subjects on the program. 
The committee would place the speakers at the proper 
tables. Members planning on attending the meeting could 
select their tables in advance. This should bring the mem- 
bers in contact with the speakers and help to break the ice 
for the industrial physicists who are comparatively young 
and timid. 

In exchange for the privilege of appearing before the 
Society the committee has a right to expect the assistance 
of the speakers in socializing the meetings. The academic 
group should be interested in cooperating, because if there 
is to be money for pure research there must be an active 
well satisfied group of industrial physicists working with 
them in an endeavor to make the results of pure research 
useful to society as a whole. 





Revise Program Material 
and Its Presentation. 
Safeguard Professional 
Interests 


. . . The lack of interest 
in American Physical Society 
affairs among indistrial phys- 
icists is evidence that some- 
thing is wrong .. . papers, 
both at the meetings and in publications, seem to be 
intended mainly for the experts and the average industrial 
physicist, who must be interested in many fields of physics, 
cannot be an expert in all. It seems that in all papers, 
even those that are largely mathematical, the results, if 
not the method, could be put in plain English. 

In general, a good paper cannot be given to a hetero- 
geneous audience in 10 minutes, and certainly all papers 
are not of equal importance. Any change would require 
that someone judge the relative importance of various 
papers... . 

Modern physics has become largely a study of electrons 
and nuclei with less attention given to anything as large 
as a molecule. At least these subjects comprise most of 
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the papers on recent programs. There are many other fields 
important to industry, for example, molecular physics and, 
to mention specific cases, molecular motions resulting in 
dipole loss and mechanical hysteresis. . . . Apparently 
many authors feel that certain types of physics papers will 
reach more interested people at other meetings. In subject 
matter, sectional American Physical Society programs 
seem to differ from those of national meetings and this 
may be significant. 

Industrial physicists are just beginning to be generally 
important and after the war there will be many of them. 
There are organizations looking out for the welfare of other 
groups such as those representing engineers, chemists, and 
the medical profession in various ways. For example, 
salary surveys are made so that an engineer or chemist can 
have some idea whether he is in the proper pay group, and 
a watchful eye is kept on legislation which may effect 
them, such as for example, the proposed Kilgore Bill. 
These particular things may not be desirable, but the 
general program is, and no existing organization is doing 
these things for physicists. 





Consider Industrialist’s I appreciate very much 
Interest in Program being given an opportunity 
to discuss the relationship of 
the American Physical Society to industrial physicists. I 
believe the chief difficulty in the situation is that the 
education of the physicist is not in accord with that de- 
manded by the industrial laboratory, and this condition 
is the prime cause of the present situation. Physicists in 
the past have been trained primarily for teaching and 
institutional research. At the present time, and for a certain 
period in the future at least, industrial institutions and 
development research laboratories can absorb a large pro- 
portion of the physics graduates of this country. 

More data are needed to indicate just what work the 
physicists of this country are now engaged in. The Amer- 
ican Institute of Physics, as you are well aware, has com- 
piled the informaticn on the present manpower situation. 
I am not aware that the Physical Society has taken any 
steps to attempt to plan programs along the lines of the 
activities of the industrial members. I suspect its plan 
has been to “‘let nature take its course.”’ I would suggest 
that the American Institute of Physics go further with its 
compilation of facts regarding its members and regarding 
the graduates in physics. I believe it will be found that 
the Physical Society now attracts to membership a smaller 
proportion of recent graduates in physics than formerly, 
say in 1910. If that is the case, you should take serious 
thought. 

Papers from industrial groups are much harder to get 
than papers from pure research laboratories. In spite of 
the fact that some very valuable scientific work comes out 
of these industrial laboratories, the situation is that these 
men are more interested in the promotion of their activities 
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along industrial lines than they are in publication. Due to 
the present attitude in Washington towards patents and 
their present tendency to invalidate patents, it will be 
more and more difficult to secure papers from industrial 
laboratories and very few papers would ever be volun- 
teered. A program of secrecy must from necessity be fol- 
lowed. 

I have yet to attend a symposium which attempted to 
tie-in the industrialist and the pure scientist which I 
thought in any way adequate. The men simply do not get 
together. It is hardly possible for the industrialist to call 
to attention his vital problems, because the things he is 
most interested in and which are most alive are no doubt 
being worked on somewhat secretly by his own company 
or group. He mostly talks about the things which have 
already been done so that the scientist is left with the 
feeling that much progress has been made and that there 
is little that he can be expected to contribute to the art, 
yet the industrialist knows that is not the case. The educa- 
tional institution is, in many cases, too much out of touch 
with the exact problems facing the industrialist so their 
attempts to cover the theoretical side of the subject also 
fall short. 





Committees on 
(1) Governmental Rela- 
tions and (2) Training 


I suggest: 

(1) The American Physical 
Society appoint a committee 
to prepare some form of ac- 
tion, possibly a session with invited speakers who will 
talk on certain phases of our Government which are of 
vital interest to the industrial physicists and their future 
such as: 


(a) The present attitude to our patent laws and the 
attitude of our central government towards the patent 
system, which as you know is an attitude tending towards 
the destruction of the value of patents. The future of the 
industrial physicist is definitely tied up with a method of 
protecting his efforts so that he and his company may 
have compensation. If this protection should cease, there 
will be no point in training men in applied physics and 
engineering. I have seen nothing about any activities of 
our society which touch this subject at all. Yet, I believe 
this subject is of great interest to any professional physicist. 


(b) The relationships of governmental laboratories and 
industrial laboratories. There is a bill, at the present time, 
in Washington which proposes to coordinate all types of 
research laboratories under governmental control. This 
you no doubt know of and I dare say that 95 percent of 
the American Physical Society, if they know of it at all, 
would be vitally opposed to such a bill, inasmuch as it 
definitely would remove the incentive motive from research. 
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(2) My second suggestion is that a committee be ap- 
pointed to investigate and recommend courses of training 
planned for the development of industrial researchers and 
substituted for the present Ph.D. course in physics. This 
committee should obtain data regarding the fields of 
physics now being worked upon by industrial physicists 
and more definite information as to the work these men 
are doing to serve as a guide for this planning. I would 
recommend to the committee that they seriously consider 
the substitution of an additional degree, call it ‘‘Doctor of 
Applied Physics” or ‘Doctor of Physical Engineering,” 
and the abolishment of the present type of Ph.D. thesis, 
substituting for that much more thorough training in 
experimental technique along the lines indicated by the 
summary of the report on the type of work that the physi- 
cists are now called upon to do. I advise this because I 
believe that very often the present experimental thesis 
does not develop research ability or experimental tech- 
nique. 

i believe that the installation of such courses and work 
will lead to publications and papers along the lines which 
tend to hold and attract industrial physicists to the Amer- 
ican Physical Society. I think these changes in our educa- 
tional program are sure to come and that the American 
Physical Society is the proper agency for the promotion of 
these or better ideas. 





Let it be understood that 
in going to the Divisional 
plan we are moving toward 


Experiment with 
Divisions on a 
Temporary Basis 
the formation of small groups 
which would work seriously through most of the duration 
of a meeting, adjourning only for the purpose of joining 
others in hearing addresses of general interest, arranged by 
the Secretary for the benefit of the whole membership. 
As soon as we do this, we gain several advantages. (1) We 
provide adequate time for real interchange of information 
between people with common interests. (2) We get away 
from the confusion of big meetings in big lecture halls. 
(3) We undoubtedly would get a much larger number of 
contributions because of the accelerating influence of Divi- 
sional Officers. Many a man would offer a paper on the 
current status of his work if he were assured of a small and 
friendly audience well acquainted with his field, when, on 
the other hand, he would hesitate to face a large crowd 
with the practical certainty that not one hearer in ten 
would have the slightest interest in what he was saying. 
Let us replace the artificial and unscientific meetings we 
have now by real working groups, held together by the 
strong bond of common interest. Only in this way can we 
recapture the delightful, unhurried atmosphere of earlier 
days. 
We should establish Divisions in several fields of 
physics at once, with the clear understanding that the 
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life of each one would be two or three years, after which it 
would continue, if successful, or would be abandoned, if 
unsuccessful, with no prejudice against its reestablishment 
at a later time. 

It is also important to realize that a very great body of 
academic physicists have turned to engineering physics 
during this conflict. It is a matter of common remark 
among these men that when they return to their regular 
places of employment they will go back with new interests 
and that many of them will continue to work on the prac- 
tical problems of industry and of government. I do not 
believe this trend will weaken the field of pure physics; I 
believe that both pure and applied physics will be bigger 
and better after the war. Let us do what we can now to 
help the practical physicists, since most of us are practical 
physicists now; and let us do what we can to have a 
greater and stronger Society, commanding the firm alle- 
giance of all physicists in the time of peace. 





Greater Promotion, 
Less Patronization 


I wish that you had ad- 
dressed me as a “physicist in 
industry” on the subject ‘‘In- 
creasing the share which physicists in industry take, etc.”’ 
instead of as an industrial physicist. 

I believe that the Society should not stress too often 
that physics is the parent science, and not attempt to 
reserve “physicist” as a term of any more distinction than 
“chemist,’’ because the attempt will more likely fail. If the 
term “professional physicist”’ is to be reserved for research 
workers only, then we can look for trouble unless and until 
colleges begin awarding degrees in engineering or applied 
physics. 

The Society will have to face the unpleasant (to some) 
fact that the use of the title Dr. as a label of distinction 
for physicists is a hindrance to full cooperation of the 
Society with engineering societies where many engineers 
of comparable training and experience do not have a 
special title of address. . . . Suppose we have in a labora- 
tory group an able electronic engineer with an E.E. and 
some years of experience, with a group of assistants in- 
cluding a physicist with a doctor’s degree and no experi- 
ence. The test of personalities is apt to be rather severe. 

To earn more money, it will be necessary to improve 
the value of the editorial matter in the Journal and in 
Review of Scientific Instruments, build up the circulation 
and advertising, hold conventions at which manufacturers 
pay for exhibit space, possibly even support a laboratory 
or two for a purpose of building up royalties on inventions. 

To a large extent, the officers of the Society are elected 
for their accomplishments in research. The marked success 
of some of the engineering societies is due to the organizing 
ability and personality of a paid manager, sometimes called 
that and sometimes secretary. The development of more 
successful schemes for promoting local chapters, divorced 
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from schools, frequently visited by one or another of the 
national officers, and thoroughly instructed and guided in 
their local management might help marvelously in the 
growth of the Society, besides helping in promoting con- 
tacts with physicists in industry. 

The book [Physics in Industry] is very convincing to 
physicists, who did not need to be convinced if they have 
been keeping themselves informed, but it must be irritating 
to engineers, at many spots, to be told, for example, that 
physicists like to thoroughly understand the fundamentals, 
with the implication that engineers do not. 





Greater Stress on 
Borderline Fields 


Your letter regarding the 
share which industrial physi- 
cists take or should take in 
the activity of our Society affairs is of great interest to me. 

There is a group here which is known as the “Physics 
Department.’’ What we do would best be called funda- 
mental work in chemical physics using infra-red and 
ultraviolet spectroscopy, x-ray spectroscopy, electron dif- 
fraction, electron microscopy, molecular ray technique, 
mass spectrography, precision calorimetry, radioactivity, 
all of which are, I should say, truly physical methods. We 
work in the fields of catalysis, reaction kinetics, lubrication, 
on analytical problems, etc. Are we chemists or are we 
physicists? Well—most of us are physicists by education 
and still like to be called physicists, although we are 
definitely working in the borderline field of chemical 
physics. We also like to work in a “Physics Department”’; 
in fact, we are proud to be physicists and we feel that 
being a physicist should not bar us from a career in the 
oil industry or the chemical industry. . . . Physicists 
actually take an increasing share of responsibilities in 
industries which, off hand, do not appear to be physical 
industries such as Bell Telephone, Signal and Switch, 
G.E., and Westinghouse. If the chemical physicists’ in- 
terest in the Physical Society (just to mention this border- 
line case) is to be preserved, the borderline must be stressed 
a good deal more than this has been done in the past. 

I should like to suggest that the Physical Society take 
cognizance of this situation by encouraging symposia and 
papers in applied physics and borderline fields. The initia- 
tive should come from the officers of the Society, and I am 
sure that the applied physicists will gladly assume the 
laborious task of arranging the symposia, etc. I personally 
am strongly opposed to a division of applied physics 
within the Physical Society, but I would favor subject 
divisions such as the proposed division of Electron Optics. 
But most of all—and I feel very strongly on this point— 
the applied or industrial physicist should not have to be 
content or should not be content to publish a piece of 
apparatus in an instrument journal or a piece of spectral 
analysis in the Review both entirely detached from the 
important purpose which his branch of physics may play 
in his industry. We should strive for a more general pre- 
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sentation of the industrial physicist and his work within 
his industry in such a way that the presentation would 
appeal to physicists as well as to scientists or engineers on 
the other side of the borderline, and possibly to company 
research executives who in this way may learn what 
physicists could contribute or are contributing to their 
industry. Well-organized symposia, properly announced 
throughout the scientific literature, will undoubtedly be the 
best medium to achieve these ends. 

Failure to see or appreciate the industrial viewpoint is 
well illustrated by the assumption that “The natural way 
is to divide physics into subjects, as for instance acoustics, 
thermionics, line spectra, etc.” This I agree is the natural 
way, from the standpoint of teaching physics and also 
from the standpoint of learning physics or cataloguing 
physics as in Science Abstracts. However, it is mot the 
only natural way of dividing the subject from other stand- 
points, including the standpoint of using physics. 

The person using physics to improve a material or a 
process must bring parts of all these fields to bear on his 
particular problem, and he likes to discuss his results and 
methods in a group whose interests and viewpoints do not 
all form an orthogonal system to his own. The need for 
this close linking of the different “natural” divisions of 
physics and their applications to specific problems is well 
illustrated in the field of metallurgy, where physics and 
chemistry and engineering are all tied together in one 
group. 





More Time, Space and 
Activity for Physicists 
in Industry 


Regarding the best solu- 
tion of this problem, I have 
no hard and fast convictions 
about what positive action 
should be taken. I do feel that merely not rejecting any 
articles submitted by industrial physicists is not enough 
and that holding symposia of invited papers is not enough. 
I think the solution will involve giving the industrial 
physicists more of a share in all of the activities of the 
Society (not only the meetings), and it will involve a 
recognition of the fact that publications by industrial 
physicists usually are the result of accumulating much 
data that cries out for publication rather than feeling that 
at least one publication a year is required for holding a 
job. One result of this is that ten minutes is not enough 
time for the presentation of any large amount of data, 
and they feel that it is not worthwhile to prepare such a 
paper. More flexibility in the time allowed for presenting 
research reports would help this situation. 

I feel that the problem will not be solved until there 
are at least as many reports of physical research published 
in the Journal of Applied Physics, as there are published, 
for example, in Industrial and Engineering Chemistry. It is 
the publication of The Journal of Chemical Physics and 
the Journal of Applied Physics which today makes the 
American Physical Society worthwhile for many physicists, 
and the council and editors are to be congratulated on 
publishing these excellent magazines. 





A National Headquarters Building for Physics 


0” July 24th the Governing Board of the American Institute of Physics approved a resolution 
authorizing the Institute to organize an appeal for funds for the purchase of a building to 
be devoted to the interests of American Physics and Physicists. An exceptionally desirable property 
at 57 East 55th Street in Manhattan had already been located and a short-term option at a favorable 
price was secured. A Physics Building Fund Committee was appointed which promptly began a 
fund raising campaign by distributing letters and other material‘ presenting the project in con- 
siderable detail to physicists and their friends. Attention was called to the desirability of acquiring 
a headquarters building as the first step in a comprehensive program for the development for physics 
outlined by the War Policy Committee of the Institute. 


The initial response to this presentation was quick and enthusiastic and by August 20th cash and 
pledges amounting to $31,816.75 had been received. On the strength of this assurance a contract 


was signed to purchase the property calling for a substantial down payment on October Ist, the 
balance of the purchase price to be paid in regular instalments during the following six months. 


To meet these contractual terms a continuing flow of subscriptions must be maintained until the 


goal of $75,000 is reached. 


Physicists and their friends who have not had an opportunity to give thoughtful consideration to 
the merits of this progressive action are urged to do so. Complete information may be obtained from 
the American Institute of Physics, 175 Fifth Avenue, New York 10, New York. 


1 Rev. Sci. Inst. 14, 241-247 (Aug. 1943). 
2 See also pages xi and xii of this issue. 
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Some Physical Procedures Used in Ceramic Research* 


By G. A. BOLE AND G. A. LOOMIS 
Engineering Experiment Station, The Ohio State University, Columbus, Ohio 


UCH engineering research procedure and 

its development into routine tests is an 
adaptation of fundamental physics and most of 
the research tools employed therein were 
originally conceived and developed by the 
physicist. This statement is particularly true of 
ceramic research, and especially so of that being 
carried on in the Engineering Experiment Station 
of The Ohio State University. The present paper 
covers some of these physical procedures as 
employed in the researches which have been or 
are now being carried on in the Station. Some of 
the physical procedures as adapted to ceramic 
research are not new while others are applications 
of quite recently developed physical apparatus 
and techniques. In describing these various tech- 
niques and their applications the detailed pro- 
cedures need not be given herein since it may be 
assumed that the physicist is familiar with them. 


DETERMINATION OF THERMAL EXPANSION 


Volume changes in ceramic products, both 
clay products and glasses, during the heating and 
cooling cycle determine largely their resistance 
to rupture when heated. A glass baking dish 
must resist thermal shock; a glaze must fit the 
body which it covers. If the glaze is in high 
tension it crazes, if in high compression it shivers 
or dunts. 

Fired ceramic systems usually consist of a 
crystalline and a glassy phase. The crystalline 
phase may be of one or more components. The 
characteristics of these phases so far as their 
reaction to thermal changes is concerned are 
recorded in the coefficient of expansion curve. 
These curves constitute significant information 
regarding the action of ceramic products over a 
given heat range. 








* Presented at the Fall Meeting of the Ohio Section, 
the American Physical Society, Mendenhall Laboratory, 
November 14, 1942. 
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The coefficient of thermal expansion is usually 
determined in one of two ways: (1) in a dila- 
tometer or (2) in an interferometer. 

A diagram of the dilatometer is shown in Fig. 1. 
The test specimen for this apparatus is a rod of 
the material under investigation about 3 inch in 
diameter and 8 inches long. This is heated in a 
fused quartz tube in an electrical resistance 
furnace at a rate of about 150°C per hour. The 
expansion of the specimen is measured by means 
of an Ames dial attached as shown to the top of 





Fic. 1. Diagram of dilatometer apparatus 
for determining thermal expansion. 
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Fic. 2. Diagram of interferometer for determining thermal expansion. 


the quartz tube and contacting the specimen 
indirectly with another smaller quartz tube or 
rod resting on top of the test specimen. The fused 
quartz has a low straight-line expansion curve 
over the area of measurement and a correction is 
made for it in the calculation of the expansion 
coefficients of the test material. 

Figure 2 shows a diagram of the interferometer. 
The test pieces for the interferometer consist of 
three tetrahedra or pyramids of the test material 
which are cut and ground to the same height 
(about } inch) or within 0.001 inch of each other. 
These are set up in the furnace as shown between 
two quartz plates which comprise the inter- 
ferometer. The light source in the viewing 
telescope is a helium-filled bulb which gives a 
monochromatic beam: When the furnace is 
heated with the test pieces and quartz plates in 
position interference fringes are observed passing 
across the field. These are counted as they pass 
a given reference point and this number with the 
rate of temperature increase (preferably 2.5°C 
per minute) provide the data for calculating the 
coefficients of expansion. The advantage of the 
interferometer over the dilatometer is the ex- 
treme accuracy of the former and the fact that 
small specimens can be used. 

Some typical thermal expansion curves are 
shown in Figs. 3 and 4. 
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DETERMINATION OF GLAZE FIT 


There are several devices for the quantitative 
measurement of glaze fit which give a direct 
indication of the relative difference in thermal 
expansion between the body and glaze. One of 
these consists of heating in an electric resistance 
furnace a two-pronged bar of the fired ceramic 
body which has been glazed on the outside only 
and measuring the spread or contraction at the 
ends of the prongs of this tuning-fork-shaped test 
piece by means of a cathetometer. 

Another method of measuring glaze fit em- 
ployed at the Station utilizes a ring of the fired 
body about 23 inches outside diameter, 2 inches 
inside diameter, and ? inch high, glazed on the 
outside face only. Reference marks are placed on 
one edge of the ring about } inch apart. Their 
distance apart is accurately measured with 
microscope fitted with a micrometer eyepiece. 
A saw cut through the ring between the reference 
marks is then made with a circular metal saw 
impregnated with diamond chips on its periph- 
ery. The distance between the reference marks 
is again measured. 
this distance 
the glaze 
the ring 


An increase or decrease in 
indicates, respectively, whether 
was under tension or compression before 
was cut. The amount of this change 
indicates the degree of tension or compression. 
It is customary in making this test by either 
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method to determine the inherent strain, if any, 
in the fired body by testing an unglazed specimen. 


THE SONIC PROPERTIES OF STRUCTURAL 
CLAY TILE 


Structural engineers desire to know the charac- 
teristics of their building units as they purchase 
these units and other building materials under 
rigid specification. Specifications and test meth- 
ods are the cornerstone on which engineering 
service operates. Most of these test procedures 
are methods and techniques originally devised 
by the physicist. 

Very little is known regarding the charac- 
teristics of ceramic building units. It is known, 
however, that these units fail if improperly 
fabricated. Failure may have its origin in the 
raw materials used, in the mechanical processing, 
the drying, or the firing of the units. The struc- 
tural materials division of the U. S. Bureau of 
Standards has been working on the problem of 
devising a method whereby the failure of im- 
properly fabricated units can be predicted. They 
have worked out a program consisting of a 
schedule for running rigidly controlled freezing 
and thawing tests and have written a tentative 
specification based on these findings. Some 
ceramic structural materials are correctly evalu- 
ated by the specification, while others are not. 
The specification does not tell either how or why 
any of the materials fail. 

When ceramic units are exposed to freezing 
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Fic. 3. Linear thermal expansion of dinnerware 
body and glaze. 
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and thawing conditions, and in fact to any tem- 
perature change, there is doubtless an internal 
friction set up. This friction may be great enough 
in time to cause strains that will ultimately 
rupture the unit. It would be desirable to know 
the trend toward failure before it has actually 
occurred. During the early stages of the changes 
which will ultimately cause rupture it would be 
desirable to be able to measure a property which 
would record this fact and especially so if it 
would record it in quantitative terms. 

A specific problem of this kind which has been 
undertaken at the Engineering Experiment Sta- 
tion of the Ohio State University is the prediction 
of failure of glazed brick and tile. Again, adap- 
tations of research tools and techniques of the 
physicist are being employed in the determina- 
tion of the sonic properties of these units. 

It is a common practice among plant men to 
tap a product to determine whether it produces 
a good “‘ring.’”’ With certain products, such as 
grinding wheels, the sonic properties are being 
measured electrically to determine their fitness 
for service. In using the sonic method to evaluate 
the physical properties of an elastic body some 
constant or value which is indicative of the 
properties of the body material must be obtained. 
One which is usually obtained in sonic determina- 
tions is the modulus of elasticity on which the 
frequency of vibration depends. This may be 
readily determined for homogeneous simple 
laboratory test pieces such as straight rectan- 
gular bars which conform to the requirements of 
ratio of thickness to length (1 to 20 or more). 
However, this method cannot be used for shapes 
with a ratio of thickness to length which does not 
conform to the determined requirements. 

The ‘‘damping properties of vibration” is an 
additional property of elastic bodies which is 
equally as important as the modulus of elasticity 
and may be determined by sonic methods on 
specimens regardless of size and shape. Excluding 
external friction the damping capacity is the 
amount of work dissipated into heat by internal 
friction in a unit volume of material during a 
cycle of vibration. A method of determining the 
specific damping capacity is by measuring the 
logarithmic decrement of vibration. This is the 
natural logarithm of any two successive ampli- 
tudes of a vibration whose amplitude decreases 
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Fic. 4. Linear thermal expansion of a clay refractory. 


exponentially and can be shown to be one-half 
the specific damping capacity. 

In the investigation of glazed brick and tile it 
is planned to determine their damping capacity 
between every fifth cycle of freezing and thawing 
in the hope that the procedure will yield results 
from which curves can be drawn that will give 
information which will indicate the genesis and 
extent of failure of the unit. 

The sonic measurement of the logarithmic dec- 
rement is a measurement of the resonance curve 
using the following formula: 


w(fo—fi) 


fo 


where e= log decrement, fp = maximum resonance 
frequency of vibration, fz and f;=frequency on 
either side of the maximum resonance frequency, 
at which points the amplitude of vibration is 1/v2 
times the amplitude at maximum resonance. 

The term (fe—f1) is also known as AF and is 
the width of the resonance band at 1/v2 of 
‘maximum amplitude. 

The resonance method must be used in order 
to measure specimens of the size and shape of 
glazed structural tile and the log decrement is 
obtained as indicative of the specimens’ prop- 
erties. 

Figure 5 shows a photograph of the complete 
apparatus used for this determination. This 
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equipment includes a beat-frequency oscillator, 
suitable vibrator, microphone pick-up, cathode- 
ray oscillograph, and the necessary amplifiers. 
The test specimen is supported on a foam 
rubber cushion so that it may be sustained in a 
free state despite the fact that its nodes of 
vibration cannot be calculated. Contact is made 
between the center of the specimen and the 
vibrator which is vibrated at the various con- 
trolled frequencies of the oscillator, and imposes 
this vibration in the specimen. The sound emitted 
by the vibrator and the vibrating specimen is 
picked up by the microphone and projected on 
the oscillograph screen. When the specimen is 
vibrated at its natural or resonant frequency a 
maximum amplitude is obtained. This point of 
maximum amplitude is shown on the screen and 
fi and f2 are obtained by changing the frequency 
to obtain the respective amplitudes at 1/v2 of 
the maximum. With these values the logarithmic 
decrement is calculated according to the formula. 


THE SONIC PROPERTIES OF GLASS 


Another study involving the determination of 
sonic properties has been made at the Engineer- 
ing Experiment Station at the request of a glass 
company. In this case it was desired to learn the 
effect on decrement of tone of variations in the 
composition of the glass, and also the effect on 
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this property of annealing and disannealing the 
glass, and that produced by altering the dimen- 
sions. 

It is customary to evaluate the quality of 
glassware such as a goblet by setting it in vibra- 
tion and observing the tone produced. A goblet 
of good quality emits a tone which is pleasing to 
the ear, principally because the vibrations are 
sustained. If the vibrations die away rapidly the 
resultant tone is unmelodic and it is concluded 
that the glass is of inferior quality. 

In the present investigation a method has 
been developed for making a quantitative meas- 
urement of the degree to which vibrations are 
sustained. This measurement was made in terms 
of the modulus of decay which is the time (in 
seconds) necessary for the amplitude or strength 
of vibration to decrease to 36.8 percent of its 
initial value (36.8 percent is the reciprocal of the 
base of the natural logarithms e, expressed in 
percent). 

The four compositions of the glasses studied 
varied mainly in their lead and soda contents 
(Table I). The work was confined to test speci- 
mens in the shape of bars } by 1 inch in cross 
section and of various lengths. The glasses were 
prepared, the test pieces pressed, annealed, and 
disannealed at the glass plant. 


TABLE I. Chemical analyses and densities of glasses used. 














Type I Type II Type III Type 1V 
SiO, 73.4 69.2 63.0 56.0 
RO; 1.0 0.1 0.1 0.3 
CaO 4.8 7.9 5.4 0.4 
MgO 3.4 2.9 1.9 0.3 
BaO 0.5 
Na2O 14.2 13.5 8.0 
KO 27 5.7 14.6 
PbO 6.3 15.7 28.5 
Density 2.6203 2.7795 3.0627 


2.4601 








A correlation of tone quality and chemical 
composition was obtained by comparing the 
moduli of decay of bars of the same size but of 
different composition. To obtain some idea of 
how the shape of a goblet affects its tone quality 
the moduli of decay of bars of different lengths 
but of the same composition were compared. 

The bar was set in vibration by a single blow 
delivered by a device actuated by a solenoid. The 
tone was picked up by a microphone and the 
pattern recorded on the screen of an oscillograph. 
A camera recorded the pattern. 

The investigation used the following equation : 
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Fic. 5. Photograph of sonic equipment used for the determination of log 
decrement. 
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where r = modulus of decay, f= natural frequency 
of vibration of bar, A,=amplitude of vibration 
at any time ¢, A,,;,=amplitude of vibration r 
cycles later. 

This equation makes it possible to determine 
the modulus of decay by measuring the ratio of 
two amplitudes of vibration, the number of 
cycles elapsed between the two amplitudes, and 
the natural frequency of vibration of the bar. 

The investigator on this study arrived at the 
following conclusions: 

(1) Within the range of lengths tested, the 
modulus of decay is directly proportional to the 
length for bars of the same cross section and 
same chemical composition. This conclusion is 
significant because it indicates the possibility of 
improving the tone quality of goblets by altering 
their design. 

(2) For bars of the same length and cross 
section but of different composition, the high 
lead content glasses have the best tone qualities. 

(3) The slopes of the lines of the graph in Fig. 
6 are such that for a given increase in length, the 
tone quality of the high lead content glasses is 
improved more than those of the low lead 
content. 
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Fic. 6. Modulus of decay of glass samples. 


(4) The effect of disannealing is to lower the 
pitch of the ring emitted by the vibrating bar 
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and to decrease the length of time the bar con- 
tinues to ring after it has been set in vibration. 


THE ELECTRON MICROSCOPE 


The electron microscope (shown diagram- 
matically in Fig. 7) is a new research tool which 
has its applications in all branches of the physical 
and biological sciences. It permits direct observa- 
tions at magnifications far beyond the range of 
the optical microscope by focusing beams of 
electrons with magnetic lenses in much the same 
manner as is done with light beams in an optical 
microscope. 

One application of this instrument is the study 
of surfaces by means of transparent replicas of 
these surfaces. Any depression in the surface is 
apparent as a dark area in a ‘‘negative’’ replica 
or as a light area in a ‘‘positive’’ replica of this 
surface, when viewed in this instrument. 

This replica method in conjunction with the 
electron microscope is being tried in connection 
with a study of the scratch hardness of glazes. 
A scratch is produced by an instrument called a 
microcharacter. The test specimen on a movable 
stage is drawn under an accurately milled dia- 
mond point. The load on this diamond point is 
varied depending on the hardness of the glaze. 
The breadth of the scratch is measured under a 
microscope. With the optical microscope the 
magnifications of the width of a scratch produced 
by the microcharacter apparatus are insufficient 
to differentiate properly between glazes which 
have a moderate degree of resistance to scratch- 
ing. Also, in the case of transparent glazes, 
erroneous measurements may be obtained be- 
cause of the uncertainty of properly focusing on 
the surface, with the high degree of magnifica- 
tion required. These scratches are only 3 to 6 
microns in width unless a sufficient weight is 
imposed on the diamond point to cause chipping 
of the glaze. The accuracy of the procedure 
depends on a track devoid of chipping. Using the 
electron microscope with magnifications up to 
50,000 diameters by photography such scratches 
may be measured in inches instead of microns. 

There are several ways of producing replicas 
of surfaces. One is to produce first a “‘negative”’ 
of the surface by condensing on the surface 
vaporized metallic silver evaporated from a 
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heating element in a vacuum chamber. This 
silver film is stripped from the surface mechan- 
ically and a “positive’’ replica produced by 
spreading a dilute solution of nitrocellulose over 
the surface, then removing the silver by dilute 
nitric acid. With this procedure difficulty is 
encountered in completely removing the silver 
without destroying the nitrocellulose film. 

Another and much easier method of producing 
a negative replica which is used for observations 
directly is to produce the film on the original 
surface with a 0.5 percent solution of polyvinyl 
formal in dioxane. This forms a tough film which 
is easily stripped off the surface without de- 
formation. 

The only apparent objection to the use of the 
electron microscope for the measurement of 
scratch hardness is that a considerable number 
of replicas will need to be made of the surface 
along each scratch as each replica must be only 
about 2 mm in diameter when mounted for 
viewing in the microscope. 


THE CYCLOTRON 


Recently the problem of measuring the thick- 
ness of the gold band on dinnerware was pre- 
sented for an immediate solution as the Federal 
government had specified that the band should 
not exceed a thickness of 1/10,000 of an inch. 
This thickness was of so small a magnitude that 
it could not be measured by any of the instru- 
ments available. A solution was therefore sought 
from the physicists and chemists. Because the 
time allowed for the solution of this problem did 
not permit of experimentation with unusual 
procedures, the technique of the chemist was 
employed, while recognizing that the procedure 
could not be expected to give an accurate result, 
under the circumstances. In this procedure, the 
gold band was dissolved off the ware and the 
gold deposited electrolytically. From the weight 
obtained the thickness of the band was calculated 
as a certain thickness, plus or minus 20 percent. 
This inaccuracy is readily explained by the facts 
that the metal deposit on the electrode was less 
than 2 milligrams, the gold lacquers used are not 
pure gold but a mixture of rare earth metals and 
the amount of these metals present is so small 
that the uncertainty of any results obtained is 
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increased. In addition to these difficulties, the 
problems of voltages and current densities re- 
quired for the depositions of the various rare 
earth metals were encountered. 




























































































Fic. 7. Diagrammatic sketch of magnetic-type electron 
microscope compared with optical microscope placed 
upside down (left). 


The cyclotron was suggested by M. L. Pool* 
as a possible solution of the problem. His pro- 
posed procedure was so simple and direct that 
it is worth recording here. He proposed to 
irradiate a section of the band with neutrons thus 
inducing an artificial radioactive gold. Simul- 
taneously a piece of pure gold leaf. would be 
bombarded with the neutrons. The radiation of 
the beta-particles from the band and the gold 
leaf would be counted in the Geiger counter. 


* Associate Professor, Department of Physics, The Ohio 
State University. 
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What the masking effect of the other components 
of the system might have been was questioned 
at the time. It is regretted that the time allowed 
for this investigation did not permit investigating 
the use of the cyclotron. 


ULTRASONIC VIBRATION AND CERAMICS 


In the purification of whiteware clays and in 
some types of clay slip casting it is desirable to 
disperse the clay in water to its ultimate particle 
size, or at least to a finer particle size than is 
possible with electrolyte alone. The magneto- 
striction oscillator is an apparatus for producing 
ultrasonic vibrations which has been found, in 
laboratory trials, to produce this higher degree 
of dispersion in the case of the more easily dis- 
persed whiteware clays. The magnetostriction 
oscillator utilizes the principle that nickel con- 
tracts when it is placed in a magnetic field and 
when it is subjected to a high frequency alter- 
nating magnetic field it produces ultrasonic 
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vibrations. The small nickel tube of the apparatus 
is inserted in the suspension to be dispersed. 

The experiments with clay slips showed that 
some of these, first dispersed with electrolyte, 
were further dispersed to a finer particle size 
when subjected to the ultrasonic vibrations. On 
the other hand, the vibration treatment alone 
without the electrolyte did not produce any finer 
particle size in most cases than did the electrolyte 
alone. A combination of the electrolyte and the 
vibration treatment produced a more dispersed 
system than did either one alone. 

Another apparatus which produces higher fre- 
quency ultrasonic vibrations than the nickel tube 
apparatus utilizes a vibrating piezo crystal 
(quartz). A million cycles per second are possible 
with it. Its capacity is less, even, than the nickel 
tube apparatus but it should have a greater dis- 
persing effect on a clay suspension. It is planned 
to try this newer technique on some clay sus- 
pensions. 


JOURNAL OF APPLIED PHYSICS 











icS 





Powdered Iron Cores 


BY ALFRED CROSSLEY 
Chicago, Illinois 


HE importance of the use of powdered iron 

cores in electrical circuits is now being 
recognized, in view of the fact that there are 
available today cores made of many different 
ferrous materials suitable for use at frequencies 
from one to 150,000 kilocycles. 

A brief historical review of the progress made 
in this new art is necessary in order that we may 
understand the reasons for the rapid advance 
made during the last decade. 

The necessity for laminated steel in low fre- 
quency alternating-current machinery and trans- 
formers to obtain desired efficient operation was 
again demonstrated by Alexanderson in the 
United States and German engineers by the use 
of 0.001”’ thick silicon steel laminations in 10,000- 
cycle alternators and 30,000-cycle frequency 
multipliers. These devices performed a useful 
purpose in radio communications, but with the 
steady increase in the use of higher frequencies, 
the laminated steel was abandoned in favor of 
air core tuning devices until about 1920 when the 
Western Electric Company introduced iron dust 
cores (Toroids) for use in telephone filter circuits, 
said filter units having far greater efficiency than 
those which included the old type thin steel 
laminations. 

These iron dust cores to the writer’s knowledge 
were not used at frequencies in excess of 60,000 
cycles for two reasons: First, the Telephone 
Company did not, in the main, require higher 
frequency services and, second, these early iron 
dust cores had particle size and insulation which 
would not provide efficient operation at higher 
frequencies (40,000 cycles and upwards). It was 
not until about 1927 that W. J. Polydoroff in- 
troduced powdered iron cores made from 
hydrogen reduced iron which provided efficient 
operation in radio tuning devices at frequencies 
as high as 1500 ke. 

In 1930 Hans Voigt in Germany introduced 
powdered iron tuning inductances, which con- 
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sisted of powdered iron particles sprayed on a 
sheet of paraffined paper which was rolled into a 
cylinder forming the core of a radio inductance. 
Later the German I. G. Farben Industries intro- 
duced a new type of powdered iron referred to as 
iron carbonyl, which due to its fineness and low 
losses is still being used in tuning devices 
operating at frequencies as high as 40,000 kc. 

The hydrogen reduced and carbonyl iron 
powder have during the last nine years been 
used in the millions of radio receivers providing 
increased efficiency and compactness of design 
unheard of in any air core inductances: 

Having briefly outlined the stages of advance- 
ment in the art, we should now consider the 
details of construction of powdered iron cores, 
mentioning first the various grades of powdered 
iron used in present-day tuning or filter devices. 

In the low frequency operating region hydro- 
gen reduced iron or permalloy materials are used. 
The hydrogen iron is produced by blowing iron 
oxide through a hydrogen furnace, which 
provides a pure iron particle having a diameter 
of from 1 to 40 microns. This iron is jagged in 
form, and, when compressed, remains in the 
shape it was compressed. Permalloy is a nickel 
iron, or molybdenum nickel iron alloy, which is 
mechanically reduced to required particle or 
grain size. It has very high permeability at low 
magnetizing field density, but loses this feature 
as the magnetic field strength is increased. Per- 
malloy when properly processed into a toroid 
form has an effective permeability of 125, 
whereas under the best conditions the hydrogen 
reduced iron toroid will provide a permeability 
as high as 50. 

In that part of the spectrum from 300 to 
40,000 ke carbonyl iron will produce excellent 
results if the proper grade is used. There are 
three grades of carbonyl namely C, E, and TH. C 
grade (4 to 10 microns diameter) is most efficient 
in the lower part of the spectrum, E (2 to 7 
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Fic. 1. Hydrogen reduced iron powder. 
Magnification 400 x. 


microns) the middle part, and TH (1 to 5 
microns) the higher frequency section. 

Carbonyl iron powder is obtained by decom- 
posing liquid iron carbonyl into the metal 
particle, which is spherical in shape and of shell- 
like structure. Figures 1 and 2 show photo- 
micrographs of hydrogen reduced and carbonyl 
iron, respectively, whereas Fig. 3 indicates the 
shell structure of carbonyl iron particles. This 
shell structure is one factor responsible for low 
eddy current losses in carbonyl powder. 

A new type of iron has just been developed 
which is of an alloy base, having the same 
properties of carbonyl TH. This iron is referred 
to as IRN-16 and is most efficient in the band 


between 5000 and 40,000 kc. Another type of 


synthetic iron powder has been produced, which 
has an apparent permeability of 2 when used 
with the best form of cylindrical winding. This 
iron is referred to as IRN-8 and is used in the 
range 30,000 to 150,000 kc. 

‘ Magnetite (lodestone) a natural mined iron 
ore when ground to required grain size, for given 
frequency range, can be substituted for hydrogen 
and carbonyl iron powder, but due to its lower 
permeability is about 75 to 90 percent the 
efficiency of the metal powders. 

In the manufacture of powdered iron cores 
the following procedures are necessary : 


1. Selection of proper iron powder for required 
frequency range. 
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2. Insulation of each iron particle from each 
other. 

3. The insulated powder is mixed with a 
binder and inserted into molds where it is 
compressed to desired form. 

4. After molding, the cores are heat treated to 
complete polymerization or aging. 


Each iron core manufacturer employs his own 
special insulation and means for processing the 
core and considers the technique a trade secret. 
This is only natural because there is no literature 
on this subject and they have learned a particular 
method of processing by the hard way, namely, 
by cut and try methods, which has taken many 
years to perfect. 

Those who require iron cores for their par- 
ticular application should specify the following 
to the iron core manufacturer so that they may 
obtain the most efficient core. 


1. Frequency range of operation. 

2. Size of inductance in millihenries. 

3. Q or efficiency of inductance. 

4. Size and type of metal container into which 
the inductance will be placed. 



































Fic. 2. Carbonyl iron “E.’’ Magnification 500 X. 
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Iron cores are made in various sizes and 
shapes, which generally are grouped in three 
classes, namely: cylindrical, toroid, and pot 
forms. Figures 4, 5, and 6 show, respectively, 
these general forms whereas Figs. 7, 8, 9, and 10 
are the types with the winding form, which are 
commonly used in high frequency tuning circuits. 

Figure 7 is the most common form and is found 
in millions of radio receivers, while the more 
efficient multiple coil unit shown in Fig. 8 is 
found in the more expensive receivers. Figure 9 
is the standard cylindrical type, the core of 
which moves in and out of the coil and, in so 
doing, tunes the radio circuit to the desired fre- 
quency when a fixed condenser is employed. This 
latter method is referred to as permeability 
tuning and is rapidly superseding the more 
expensive variable condenser type of tuning. 
Figure 10 is the pot type wherein the coil is com- 
pletely enclosed. The latter type when used with 
proper type of litzendraht cable is the most effi- 
cient type of high frequency inductance for use 
at frequencies higher than 50 kc. The writer has 
constructed a 200-microhenry inductance of this 
type which has a Q of 420 at 1000 kc. 

Space does not permit a complete discussion of 
iron core tuning devices but it is permissible to 
mention a few outstanding applications where 
powdered iron cores produce remarkable results. 





Fic. 3. Carbonyl iron “‘E.”’ No. 43-s. 
Magnification 2000 x. 
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Fics. 4-6. 


The first device is the cylindrical aircraft 
direction finding loop antenna (Polydoroff patent 
No. 2,266,262) using powdered iron cores. With 
this construction the same pick-up voltage is 
obtained with a 23” diameter iron core loop as 
was possible with a 9” air core loop. The smaller 
iron core loop reduces considerably wind re- 
sistance or drag over that of the larger air core 
type and in addition provides a better minimum 
or zero signal for the right angle position of the 
plane of the loop to the advancing wave front. 

In this type of loop antenna the increase in 
efficiency is due principally to the increase in 
permeability of the media caused by the presence 
of powdered iron in the field of the coil. The 
reduction in size also cuts down the so-called 
“antenna effect’”’ thereby producing a better 
minimum or zero (right angle) signal. 

It is only natural to assume that all future 
aircraft loop antennas (direction finding or radio 
compass) will be of the iron core type for the 
obvious reasons previously mentioned. 

The second application is a new device which 
will revolutionize audio- and radiofrequency 
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coupling systems. This device provides means for 
varying coupling from zero to a desired maximum 
value and in addition the phase relationship of 
said inductive coupling can be adjusted to be 
aiding or opposing the electrostatic coupling 
inherent in the system being coupled. 

Figure 11A is a schematic drawing of the con- 
ventional coupling transformer used in the radio 
receivers, comprising a primary winding with its 
tuning capacitance and a secondary winding also 
with the tuning capacitance both of which permit 
transfer of power from the plate circuit of one 
vacuum tube to the grid circuit of the succeeding 
tube. It is common practice to couple mag- 
netically one coil to another by moving one coil 
to a given distance from the other so as to 
obtain critical coupling. Increasing the coupling 
beyond critical coupling, introduces complica- 
tions, because the mutual inductance starts in- 
creasing both primary and secondary inductance 
and both circuits are thereby tuned to a lower 
frequency. 

Figure 11B shows a means for compensating 
for increase in mutual inductance by decreasing 
automatically the inductance of the primary and 
secondary when the mutual inductance or coup- 
ling is increased or vice versa. In this circuit 
three powdered iron cores linked together are 
associated with the following respective parts of 
the circuit: first, a small coil in series with the 
primary; second, a coupling coil; and third, 
another small coil in the secondary circuit. The 
movement of the iron core into the field of the 
coupling coil also produces a movement out of 
the field of both primary and secondary coils 
providing an increase in coupling or mutual 
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inductance, but at the same time a decrease in 
inductance of the primary and secondary in- 
ductances thus keeping constant the total in- 
ductance of both primary and secondary circuits. 

The maintenance of a constant inductance in 
both circuits with variation in coupling between 
circuits maintains the tuning of both circuits to 
a fixed frequency and in addition provides means 
for obtaining extremely sharp tuning or any 
intermediate coupling up to the over-coupled 
condition necessary for high fidelity reception. 
With such a system the side band acceptance is 
uniform because there are no distorted side band 
characteristics. Figure 12 indicates the receptivity 
characteristics of the system for the varying 
degrees of coupling. 

Figure 11C is an improvement electrically and 
economically over arrangement shown in Fig. 
11B. In this figure one coupling coil system and 
one iron core do the work of three previously 
used. The coil consists of a Bakelite tube upon 
which is wound bifilarly four coils and inside of 
the tube is a powdered iron core having an ad- 
justing rod. The winding connections of one of 
the secondary coils are reversed thereby pro- 
viding the novel condition whereby the move- 
ment of the iron core can produce any desired 
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degree of coupling from zero to maximum with 
the additional feature of phase changing. 

Zero coupling is obtained by moving the core 
to the center of the tube whereby the negative 
mutual obtained with the reversed primary coil 
is counterbalanced by the positive mutual or 
coupling of the lower coil system. Moving the 
core upwards produces an increase in coupling 
until maximum is obtained when the core is 
completely enclosed by the upper coil system. 
Moving the core downwards from the center 
position also increases the coupling but this 
coupling is 180 degrees out of phase to the 
previously mentioned coupling. 

Throughout the entire movement of the core 
for any degree of coupling the inductance of both 
the primary and secondary remains constant. 

This coupling system can be applied to oscil- 
lating systems for positive control of feedback. 
It also can be used with amplifier and filter 
circuits where negative or positive control of 
coupling or feedback is required. 

In concluding it is permissible to say that we 
have, in the powdered iron core, a tool which 
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KILOCYCLES OFF RESONANCE 


Fic. 12. Variable coupling interstage transformer. Fre- 
quency 456 kc. 1.5 millihenries, primary and secondary. 
No. 1, under-coupled, gain 125; No. 2, under-coupled, 
gain 263; No. 3, optimum coupled, gain 286; No. 4, over- 
coupled, gain 286; No. 5, over-coupled, gain 250; No. 6, 
over-coupled, gain 112. 


permits us to obtain three important conditions: 
first, increase of inductance with increased effi- 
ciency ; second, tuning of or coupling changes in 
circuits by movement of iron core; and third, 
increase in signal voltage in loop antennas by 
increase in permeability of the media. These 
functions, when properly applied, will open up 
new fields for research in both audio- and radio- 
frequency circuits. 
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Frictional Phenomena. XV 


By ANDREW GEMANT 
The Detroit Edison Company, Detroit, Michigan 


Chapter XV. External Friction of Solids 


Abstract 


N this chapter external friction between sliding 

solids is discussed. The fundamental mecha- 
nism, which is based on the geometric irregularity 
of macroscopically plane surfaces, is discussed 
first. A review of the experimental facts on the 
dependence of the friction coefficient on pressure, 
speed, and temperature reveals the plastic 
deformation of the surface irregularities as a 
complicating factor during sliding. Further ex- 
perimental material on the so-called slip-stick 
motion, and on observations of traces of soft 
metals left upon harder metals sliding against 
them indicate abrasion as another complicating 
factor in external friction. 


57. Basic Mechanism of External Friction 


The last field to be dealt with in this mono- 
graph is the friction experienced when two solid 
bodies are shifted against each other. This type 
is certainly the first that came to any observers’ 
attention as it is the most obvious and most 
frequent in both everyday and scientific proc- 
esses. Other fields, such as liquid viscosity, were 
conceived and explored at a later time, and were 
distinguished as internal friction, in contradis- 
tinction to the external friction that occurs at 
the boundary between two solids in relative 
motion. 

In spite of this primary nature of external 
friction it has been placed at the end of this 
treatise, the sequence in which, as has already 
been stated, is one of increasing complexity. And 
indeed, external friction, which was observed 
first, is apparently the most complicated as to 
its mechanism. While it has been possible to 
present a fairly satisfactory, although often far 
from complete, theory of gaseous, liquid, non- 
Newtonian, plastic, and internal solid friction, 
this is by no means possible in the field of external 
solid friction. There simply does not exist any 
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generally accepted theory; all we have are only 
first and very fragmentary attempts along this 
line. 

The procedure to be followed in this chapter 
will be this: First, a very elementary mechanism 
underlying external friction will be presented. 
This mechanism is most likely correct, although 
really it is only a first step to a theory since 
probably a large number of complicating proc- 
esses play a more or less important role in 
individual cases. Next, the existing experimental 
material, (which, too, by the way, is very incom- 
plete, and at the best only semi-quantitative) will 
be presented. In discussing the experimental 
facts, the various additional complications in the 
mechanism of friction will be considered in such 
a manner as to give an explanation for the various 
experimental observations. This procedure is 
about the best that can be followed with our 
present knowledge. 

The basic mechanism of external friction 
originates from Coulomb. It is known that all 
solid surfaces are rough, the roughness being 
either macroscopic (dimension of asperities or 
protuberances ~10-* cm), or microscopic (~ 107% 
to 10-* cm), or even submicroscopic (~ 10-5 cm). 
The dimensions of the asperities along a surface 
are probably randomly distributed over certain 
average values, but for sake of simplicity we shall 
operate in the following with these average di- 
mensions. Two horizontal surfaces, A and B, 
sliding on each other are shown schematically in 
Fig. 142, A moving to the right, and B being at 


Fic. 142. Schematic diagram of two bodies sliding 
upon each other. 
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rest. The weight of A, shown by the arrow P, 
presses A against B with a pressure p=P/a, if 
a is the contact area. 

If A is moved with a constant speed to the 
right, the hills of the top surface will drop peri- 
odically into the valleys of the lower surface. If 
the average distance of the asperities is \ cm, 
then the frequency of this periodic motion will 
be 1/A per cm of path covered. Because of the 
irregularity of the surface pattern, the phase of 
this periodic motion will be different for the 
various asperities, even if the period is approxi- 
mately the same for all. Such a coexistence of 
different phases of motion among the asperities 
is obviously possible only if deformations, 
whether elastic or plastic, of the asperities take 
place. As a first rough approximation, however, 
these deformations will be disregarded. 

To raise a hill out of a valley into which it has 
dropped, in order that horizontal motion may be 
maintained, will require a certain force applied 
horizontally. The magnitude of this force per 
surface unit Fy can be estimated by considering 
a single hill of A with a base area \?. The hori- 
zontal force acting on it is Fo’, the vertical 
pressure is pd*. Consider a motion of this hill as 
it is lifted from a valley of B to a hill-top. Such 
a motion will in reality not require a constant 
force, since the latter depends on the slope of the 
hill which varies from point to point. If, however, 
F)\* represents the average of the force required 
to move a hill from one valley to the next, its 
work is given by FoA*, since \ is the total path 
for each period of the motion. The potential 
energy obtained is pd*h, if h is the height of the 
hill. Hence, equating the two, we have 


Fy=(h/d)p. (166) 


Here each side of the equation has been divided 
by \*, in other words, the relation, valid for a 
single hill, has been averaged over the surface 
unit. 

The question immediately arises: What hap- 
pens with the gained energy during the downward 
motion that follows each lifting? When a hill of 
A passes over the top of a hill of surface B, it 
will be pushed down into the next valley owing 
to the pressure p. This impact is now assumed 
to be essentially inelastic, i.e., the potential 
energy ph will be converted into heat. A 
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mechanism of this kind, in which the work done 
by the external force is dissipated as heat, must 
lead to the observation of a friction force, the 
magnitude of which is given by Eq. (166). 

It is quite obvious that this crude picture of 
external friction based solely upon geometric 
consideration of the sliding surfaces, can be con- 
sidered as a framework at best. The asperities are 
certainly not rigid, as assumed here, but, as 
already mentioned, will undergo elastic deforma- 
tions, both in the direction of and normal to the 
surface. Several other complicating factors will 
become evident from the next section. The con- 
ception of external friction as a superposition of 
a number of factors upon a basic mechanism was 


also emphasized recently by Claypoole and 
Cook.! 


58. Influence of Experimental Factors on 
External Friction 


(1) Effect of Pressure and Area of Contact 


One of the fundamental laws of external 
friction, first stated by Amontons, relates the 
total frictional force F, acting between two 
bodies, to the weight P of the top body or, more 
generally, to the total pressure, acting normal to 
the boundary. The relation states that the two 
quantities are proportional and is expressed by 


F=pP. (167) 


The proportionality factor is called the coefficient 
of external friction, and it is a quantity depend- 
ing only upon the nature of the two surfaces in 
question. 

This figure has not great significance as to its 
absolute value, since it depends on many factors, 


TABLE XIX. Coefficient of external friction yu for various 
materials, after Tichvinsky, and Schnurmann. 











~ Combination of materials “ 
Carbon —Glass 0.18 
Copper —NMild steel 0.36 
Garnet —Mild steel 0.38 
Glass —Glass 0.40 
Hard steel—Hard steel 0.42 
Cadmium —Mild steel 0.46 
Ebonite —Glass 0.53 
Mild steel—Mild steel 0.57 
Copper —Copper 0.60 
Nickel —Mild steel 0.66 
Cadmium —Cadmium 0.80 
Aluminum—Aluminum 1.4 
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but for sake of illustration a few data, as 
measured by Tichvinsky? and by Schnurmann,? 
for dry surfaces are given in Table XIX. It can 
be seen that the figures are of the order 0.2-0.7, 
but figures higher than 1.0 are occasionally ob- 
served. Figure 143, also after Tichvinsky, shows 
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Fic. 143. Friction coefficient » vs. load P, for various 
materials: 1. Nickel—mild steel; 2. Hard steel—hard steel; 
3. Glass—glass; 4. Carbon—glass. (After Tichvinsky.) 


that the coefficient is, indeed, independent of the 
load, proving that Eq. (167) is correct. 

Amontons’ law contains the additional state- 
ment that the total frictional force F for a given 
total load is independent of area of contact, a 
statement often verified experimentally. 

Both these statements follow directly from Eq. 
(166) of the previous section. Multiplying both 
sides of this latter by the area of contact a, Eq. 
(167) immediately obtains, with 


p=h/n. (168) 


Thus yu, the friction coefficient, depends only 
upon the nature of the surface, as has been found 
experimentally. The relation (168) is, of course, 
just as crude as was Eq. (166). 

In practical cases the surface of contact 
between two bodies is often curved. The actual 
area of contact in such cases then increases with 
the total weight P. Yet, Eq. (167) is valid: The 
frictional force is proportional to the total 
pressure and is not affected by the variation of 
the contact area with pressure. 


(2) Velocity of Sliding 


Observations show that the coefficient of fric- 
tion for low relative speeds is approximately 
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independent of the speed. This rule was verified 
by Tichvinsky, among others. He took decelera- 
tion curves which, when analyzed, give directly 
the friction as a function of speed. If the two are 
independent, deceleration is constant. He found 
the latter to be true between relative speeds of 
60 and 600 cm per sec. 

The rule appears to break down for higher 
speeds, in that » decreases with increasing ve- 
locity. The speed limit at which the constancy 
of u ceases appears to vary considerably with the 
materials investigated. Some data are shown in 
Fig. 144, as compiled by Stanton.‘ The coef- 
ficient is shown as a function of relative speed in 
mph (1 mph=44.7 cm/sec.). Curve 1 refers to 
cast iron brake blocks on steel wagon wheels, 
curve 2 to locomotive driving wheels on steel 
rails in foggy weather. (This last qualification is 
essential, since water vapor introduces boundary 
lubrication ; see Section 59 (2), also Chapter X.) 
It should be added that curve 2 refers to so-called 
rolling friction ; although important from a prac- 
tical standpoint, this type will not be dealt with 
in this monograph. 

It is customary to differentiate between static 
and dynamic friction coefficient, the first cor- 
responding to the force just necessary to over- 
come friction and start sliding, the second to the 
force necessary to maintain motion. It appears 
that under ideal conditions (perfectly dry, hard 
surfaces) the static coefficient coincides with the 
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Fic. 144. Friction coefficient yu vs. relative speed: 1. Cast 
iron brake blocks on steel wheels; 2. Steel wheels on steel 
rails in foggy weather. (After Stanton.) 
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dynamic coefficient for low speeds. In most prac- 
tical cases, however, an abrupt drop in the value 
of u takes place, as soon as motion starts. This 
effect is probably a secondary one, due to abra- 
sion of the tips of the hills. 








Fic. 145. Angle of repose. 


If the force acting upon a body that is detained 
by external friction increases from zero, no 
motion occurs until the force corresponding to 
the static coefficient is reached. This force may 
be looked upon as a yield value, which is very 
well defined in this case, in contradistinction to 
plastic motion where it was shown to depend on 
the sensitivity of the measuring device. 

The static coefficient is related to the so-called 
angle of repose, which is obtained by tilting the 
two bodies until sliding sets in (Fig. 145). Let us 
denote this angle by a. With Po= weight of the 
top body, the total pressure acting upon the 
lower surface is Py) cosa@ and the total driving 
force is Po sin a. Substituting Po cos a for P and 
Py) sin a for F in Eq. (167), we obtain: 


p=tan a. (169) 


The angle a is called angle of repose; by means 
of Eq. (169) it allows the static value of u to be 
determined. 

In the light of the discussion above it is evident 
that accurate figures on » should be specified 
whether the static or a kinetic value is meant 
and, in the latter case, at what relative speed. 

The experimental fact that the frictional force 
is initially independent of the speed and de- 
creases at high values of the latter is in quali- 
tative agreement with the Coulomb mechanism. 
It can be seen that Eq. (166) does not show the 
speed explicitly. However, it can be shown that 
for high speeds a decrease in the effective value 
of h, and thus of yw, might be expected. For small 
speeds the asperities will drop to the bottom of 
the valleys, whereas for high speeds the time will 
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not be sufficient for a drop of this extent to take 
place and the asperities will hit instead the 
ascending branch of the next hill. Thus the 
effective depth of the valleys will be lowered. 
While this explanation seems to be correct 
qualitatively, it does not hold quantitatively. A 
decrease of » can be expected when the order of 
speed becomes 


uo~d(g/h)$ (170) 


(g= gravitational constant). With \ and h/ of the 
order 10-* cm the value of %) would be of the 
order of 1 cm/sec., whereas the experimental 
data in certain cases indicate a higher order than 
that. This also shows that other processes apart 
from a mere lifting against the normal pressure 
take place during sliding ; these will be discussed 
below. 

It should be added that there are exceptions 
to the above-mentioned typical variation of u 
with speed, particularly in the case of rubber. 
While the friction coefficient at speeds of less 
than 10-* cm/sec. is less than the static coeffi- 
cient, for higher speeds the coefficient increases 
with the speed. This behavior is brought out 
clearly in a recent paper by Roth et al.5 This 
exceptional behavior of rubber is not surprising 
in view of its low elastic modulus, making the 
simple concept of lifting of asperities of minor 
significance in this case. 

Interesting results are obtained if the de- 
pendence on speed is considered in connection 
with the analogous effects in other frictional 
processes previously discussed in this monograph. 
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Fic. 146. Angle of repose vs. temperature for cadmium— 
cadmium interface, according to Schnurmann. 
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However, it will be more useful to discuss this 
matter after the temperature effect has been 
considered. 


(3) Effect of Temperature 


Although the question of the temperature 
influence is of considerable practical and theoret- 
ical importance, the information upon this point 
is very scanty. Schnurmann® recently inves- 
tigated this question by using materials that 
were baked out at 160°C and then vacuum- 
treated in order to remove absorbed layers of 
vapors. Some of his data are shown in Fig. 146. 
They refer to a cadmium-cadmium interface, and 
show the angle of repose a as a function of 
temperature. While the angle (and yu) appear to 
be fairly constant below room temperature, 
above this up to 150°C a slow decrease takes 
place, the value of da/dt being —0.08. It should 
be recalled that cadmium has the low melting 
point of 321°C; thus the range of decreasing yp is 
not too far from the melting point. Other metals 
of higher melting point may not show a decrease 
near room remperature. 

Merchant® investigated polished pairs of like 
metals; when they were tested below the an- 
nealing range, » was found to be independent of 
temperature. 

The simple model explained above does not 
indicate any influence of temperature on external 
friction. Thus, it appears that the model might 
approach reality at lower temperatures, but has 


TABLE XX. Summary of various frictional processes. 


State of material Mechanism of friction 





to be modified at higher ones. The modification 
consists in introducing the plasticity of both or 
at least one of the sliding bodies as a factor. The 
asperities do not have a constant shape, but may 
be plastically deformed under the action of the 
tangential force. The chief result of this mecha- 
nism is a reduction of the effective value of h, the 
height of the asperities and, hence, of the friction 
coefficient pw. An accurate evaluation of this 
influence of plasticity on external friction has 
not yet been carried out. It will become effective 
if at least one of the materials has a finite plas- 
ticity—for instance, if it is nearing its melting 
point. The effect obviously increases with in- 
creasing temperature (decreasing plastic re- 
sistance) in agreement with the behavior shown 
in Fig. 146. 

At this point of the discussion it will be useful 
to review briefly the different friction mechanisms 
encountered in this monograph. Table XX gives 
a summary. The first column shows the state of 
the materials. The last item in that column con- 
tains the specification ‘‘dry’’ in order to dif- 
ferentiate it from lubricated surfaces (see Sect. 
59 (2)). The second column shows the type of 
mechanism underlying friction, the third the rela- 
tion between stress and rate of shear, and the 
fourth the relation between frictional force and 
temperature. 

The chief point to be emphasized is that as 
long as viscosity is controlled by the kinetic 
energy of molecules, the stress increases with rate 





Frictional force vs. 
temperature 


Stress vs. rate 
of shear 





Gaseous Molecular, purely kinetic 


Liquid Molecular, potential barriers 


vs. kinetic energy 


Liquid, non-Newtonian Molecular + submicroscopic, 


potential barriers vs. kinetic 


energy 


Solid, plastic Molecular + submicroscopic, 


potential barriers vs. kinetic 


energy 


Solid, surface dry Submicroscopic + micro- 


scopic, potential barriers 





Increasing, according to 3 


power 


Increasing, linearly 


Increasing, linearly Decreasing, exponentially 


Increasing, less 


Decreasing, exponentially 
than linearly 


Increasing, loga- Decreasing, exponentially 


rithmically 

Independent, Independent, secondarily 
secondarily de- decreasing 
creasing 
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of shear and, obviously, there will be a pro- 
nounced variation of viscosity with temperature. 
An explanation for the first of the above twostate- 
ments can be formulated generally in the follow- 
ing way. Kinetic energy involves the validity of 
statistical laws, hence the number of moving 
molecules (rate of shear) will vary according to 
the Boltzmann function, the potential energy in 
which is both internal and external (stress). 
When kinetic molecular energies cease to play a 
role, as in external solid friction, there will be no 
variation of stress with rate of shear (unless 
secondarily), and also no variation of the stress 
with temperature (unless secondarily). The 
process is then controlled by potential barriers 
alone and the rate of shear does not enter into 
the picture. These interesting relations are best 
illustrated by the plastic resistance of solids, 
which occupies a transitory position between 
frictional forces controlled by kinetic and poten- 
tial energies, respectively. It was shown in 
Section 42 how the plastic stress approaches a 
constant value as the elastic potential energy 
becomes commensurable with the potential 
energy of the barriers; in this stage the kinetic 
energy has lost its significance. 

The trend shown by the stress vs. rate of shear 
relation for various states of matter is shown 
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Fic. 147. Schematic diagram of stress vs. rate of shear, 
for various states of matter. 


diagrammatically in Fig. 147. The coordinate 
system has to be shifted along the curve in order 
to show the behavior of the various states of 
matter, indicated along the curve. In order to 
include yield values, in particular the static 
friction coefficient, the origin of the coordinate 
system has to be placed at the three respective 
circles shown in the diagram. 
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59. Further Experimental Observations on 
External Friction 


(1) The Springiness of the Sliding Parts 


Interesting effects are observed if both parts 
are not rigidly mounted, one of them possessing 
a certain flexibility. In order to fix our ideas let 
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Fic. 148. Coefficient of friction vs. time in a stick-slip 
process for bearing metal on steel. Motion from left to 
right. After photogram obtained by Bowden et al. 


us assume that the lower surface is rigidly at- 
tached to a mechanism (a piston, for instance) 
and has a constant speed. The upper surface, 
which is at rest, may be mounted rigidly or not. 
In the former case the relative speed is strictly 
constant at any time. In the second, when the 
upper body is mounted by means of springs, a 
discontinuity of the motion is observed. These 


are the results of certain vibrations exhibited by 


the top body. Such vibrations, which are caused 
by a steady motion controlled by external fric- 
tion, have been known for a long time and have 
practical significance. This latter aspect will be 
discussed in Chapter XVI, while some basic 
experiments, with their results, will be presented 
in this section. These were carried out by Bowden 
et al.” ® and more recently by Muskat and col- 
laborators.°® 

If the motion of the top body is registered by 
means of a moving film camera, pictures like that 
shown in Fig. 148 are obtained. This figure was 
drawn after a photogram obtained by Bowden. 
The abscissae are units of time (seconds) and the 
ordinates are displacements that have been con- 
verted into force, and, by dividing by the load, 
into friction coefficient. 

The top surface is first dragged along with the 
lower one. This can be seen from the slope of the 
line, which, in terms of displacement, is found 
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to be identical with the speed of the lower surface. 
The elastic restoring force of the top body then 
increases, obviously until the static friction force 
is reached (in the figure at about 0.5). What 
happens next depends upon the friction force vs. 
velocity characteristic. 

If the force were independent of, or were to 
increase with, the velocity (represented by a 
so-called positive characteristic), then a stable 
equilibrium would be attained. In case of inde- 
pendence, a certain deflection of the top spring 
would result independently of the relative speed ; 
in case of a positive characteristic, the deflection 
would increase with the speed, but the equilib- 
rium would be stable and the motion stationary 
in both cases. 

If, however, the friction force decreases with 
increasing velocity, as represented by a so-called 
negative characteristic (see the third section of 
the curve in Fig. 147), then no stationary equi- 
librium is possible. As a relative speed develops, 
the friction force drops, resulting in an accelera- 
tion of the top body toward its normal position. 
Owing to this impulse the latter will perform 
vibrations in its natural frequency, which will 
decay according to the total damping present. 
The ‘natural frequency was high in Bowden’s 
case ; thus these vibrations are not visible in Fig. 
148. In one figure published by Bowden, however, 
these oscillations were clearly visible. When the 
relative speed again drops to a low value, the 
friction force increases and the top body again is 
dragged along, as seen from the figure. The whole 
process is repeated again. Hence, as has been 
explained, this remarkable periodicity, termed 
stick-slip, is caused by the negative characteristic 
of force vs. velocity. 

It should be added that if the vibrational 
velocity of the elastically supported body and 
the velocity of the continuously driven body 
‘become of the same order of magnitude, the 
resulting stick-slip process will develop into 
undamped, self-excited vibrations, corresponding 
to the natural frequency of the elastic body (see 
Section 63). 

Useful information is obtained from measure- 
ments of the contact resistance between the two 
bodies. They indicate that during the stick-phase, 
as the force of friction increases, the points of 
contact undergo some changes. This is brought 
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out clearly on the two diagrams shown in Fig. 
149, drawn after photographs of Bowden. They 
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Fic. 149. Coefficient of friction and electrical conductance 
vs. time in a stick-slip process. After photograms obtained 
by Bowden et al. (a) Constantan—steel. (b) Tin—steel. 


show simultaneous variations with time of the 
friction coefficient and the electrical conductance 
of the contact. Diagram (a) refers to constantan- 
steel, (b) to tin-steel. Increasing conductance 
points to increasing contact area of the protuber- 
ances. A change in the contact area probably re- 
sults from a deformation of the points of contact. 

In the first of the two cases shown, increasing 
frictional force flattens out the points of contact 
and conductance increases. When the slip takes 
place, the conductance drops rapidly to the 
original value. 

In the second case the drop in conductance 
occurs toward the end of the stick-phase. Bowden 
explains this by assuming that welding of the 
contact points takes place, and that these minute 
welds become extended toward the end of the 
stick-period, thus causing a drop in the con- 
ductance. 

The question of welding is a controversial 
point in the theory of friction. It appears'® that 
a raise in temperature to a sufficiently high level 
for a welding to take place is quite exceptional. 
Such a temperature rise might happen if the 
sliding velocities were particularly high, or if the 
melting points of the components were par- 
ticularly low. Because of the high thermal con- 
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ductivity of metals, it is not very likely that 
locally high temperatures could be obtained. 
Nevertheless, the possibility of this complication 
should be kept in mind, for instance, at very high 
pressures at which the sliding surfaces stick. 
This so-called seizure pressure depends upon the 
finish of the surfaces.'' However, welding cannot 
be considered as a generally present feature in 
friction; it certainly does not take place during 
friction of non-metallic surfaces, like wood. 

Even if the theory of a welding of contact 
points does not hold generally, there is a possi- 
bility of contact points of either of the two bodies 
approaching the other surface to such a degree 
that molecular attraction forces between the two 
bodies become active. It is known that such 
forces extend only to a distance of the order 10-7 
cm, but an approach to this degree might take 
place during the stick-period. If this is admitted, 
then the simple Coulomb model, based on geo- 
metric considerations, has to be modified by 
taking into account the molecular attraction as 
part of the frictional force. According to this idea 
the molecular attraction forces have to be over- 
come each time a barrier is passed. Tomlinson” 
worked out a mathematical theory along these 
lines. It is not easy on the basis of available ex- 
perimental material to decide whether such a 
mechanism plays an essential role or not. Cer- 
tainly in those cases in which there is little or no 
dependence of friction on temperature this 
mechanism cannot be too important. 

Whatever the role of local welding or molecular 
attraction forces in certain cases, their presence 
cannot be concluded from the mere fact of the 
stick-slip process. The negative friction vs. speed 
characteristic is, as pointed out above, sufficient 
to produce the effect. 


(2) The Nature of the Sliding Surfaces 


Further interesting observations are made if 
the sliding surfaces are investigated microscopic- 
ally after sliding has taken place. 

If the two materials have widely differing 
hardnesses (as in the case bearing metal and 
steel), then the softer material often exhibits 
characteristic parallel lines upon its surface in 
the direction of sliding. These lines can be inter- 
preted as grooves cut by the harder material on 
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the softer. Such lines are observed not only with 
two dissimilar metals, but also with glass-metal 
combinations, as shown by Schnurmann.’® The 
harder of the two metals shows characteristic 
traces of the softer metal upon its surface after 
sliding. 

These observations indicate the most radical 
departure from the simple geometric model. They 
clearly show that asperities, when lifted upon 
each other, do not always glide over the top, nor 
is the only possible complication a deformation 
but that a breaking away of the protuberances 
often also takes place. The work spent against 
the frictional force in such a case is used partly 
in the creation of new surfaces, this being to a 
certain extent an exception to the general rule 
that work against friction forces always appears 
as heat. 

When a harder metal slides upon a softer one, 
it is chiefly the asperities of the latter that break 
off and then adhere to the harder; hence the 
grooves appearing upon the softer metal and the 
traces appearing on the harder. The final effect 
will be that the small amount of the softer metal 
that has been abraded now acts as a sort of 
lubricant between the two surfaces, a rather 
beneficial effect from the practical standpoint. 
The greater the difference in hardness between 
the metals the more pronounced will be the 
effect. Soft metals of this type, copper, tin, or 
lead alloys, are used as bearing metals."* 

However, if the two surfaces are the same, each 
will abrade the other, no soft lubricant will be 
produced, and the wear of such a pair of equal 
metals will reach much higher grades than if one 
of the pair is softer. The practical implication of 
this conclusion is obvious. 

It is not only the hardness of the sliding sur- 
faces that affects the friction. Schnurmann" 
stresses the fact that an adsorbed layer of water 
or other impurity is often present upon solids and 
will considerably influence the friction coefficient. 
This effect'® was discussed in some detail in 
Section 40 of Chapter X. While it deserves con- 
sideration from a practical standpoint, it is be- 
lieved that the adsorbed layers are not generally 
characteristic for the mechanism occurring be- 
tween unlubricated, dry surfaces, which for this 
reason is often called dry friction. This latter is 
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more adequately described by the three ele- 
mentary processes taking place at the pro- 
tuberances: lifting, deformation, and abrasion. 
Other complications, such as welding or bound- 
ary lubrication by absorbed layers, are rather 





incidental. In the latter group of complications 
belongs also a component of friction due to 
contact electrification, observed by Schnurmann 
and Warlow-Davies'® on electrically insulated 
surfaces—even with a pair of metals. 
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214 (1939). 
HE calculation reported on page 229 should 
be corrected to read as follows: 
2 Percent Extension 
b=2.35/4=0.59 
0.35g/d=A —0.59 log (10*+1) 
A=2.12 
S=2.12—0.59 log (t+1) 
S=2.12—0.59 log (10-?+1) 
S=2.12g¢/d 
yg 2:12 1.33 K 10°F 
poi 0.75 
= 3.810" dynes/cm’. 
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4 Percent Extension 


b=2.25/4=0.56 

0.55¢/d =A —0.56 log (10°+1) 
A =2.23 
S= 2.23 —0.56 log (t+1) 
S= 2.23 —0.56(10?+1) 





S=2.23¢/d 
2.23X1. p 
4% 


=7.4X 10" dynes/cm*. 


* These values do not represent instantaneous moduli 
but rather an average of the composite moduli in effect 
during the respective extension. 

+ 1 gram/denier = 1.33 X 10° dynes/cm?*. 
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Résumé of Recent Research 


New Books 








The Structure of 
Liquid Mercury 


Many liquids have 
now been investigated 
by means of x-ray dif- 
fraction techniques similar to the Debye-Sherrer 
system for crystals, in attempts to determine the 
average arrangement of the atoms in the liquid. 
Monatomic liquids are particularly interesting 
since they may be treated more exactly and 
should, in general, give simpler results. 

Campbell and Hildebrand! have investigated 
the structure of liquid mercury using rocksalt 
monochromated Mo Ka radiation. Patterns were 
obtained under atmospheric pressure at — 38°, 
0°, 50°, 100°, 150°, and 200° C. The x-rays were 
diffracted by a horizontal mercury surface at an 
angle of incidence of 4°45’ and recorded on a 
quadrant of photographic film a distance of 7.00 
cm from the sample. 

Analysis of the intensity patterns showed that 
there are approximately 6 (6.00 to 5.25) atoms 
in the first coordination layer of the liquid at an 
average distance of 3.00A from the central atom. 
At the lower temperatures a marked second 
coordination layer of 5 to 6 atoms was found in 
the neighborhood of 4A. This second shell gradu- 
ally disappears as the temperature is raised. 
Further concentrations were observed at dis- 
tances of 6A and 7A from the central atom. 

The first coordination layer corresponds closely 
to the arrangement in the solid where 6 atoms 
occur at 3.00A with 6 more at 3.4A. This liquid 
coordination number of 6 is considerably lower 
than reported values for other liquefied mon- 
atomic elements which may account for some 
of the very unusual properties of liquid mercury. 

As has been predicted from the general con- 
ception of liquid structure the average distance 
to the first coordination layer varies only slightly 
with temperature and the distance of closest 
approach of the atoms decreases with rise in 
temperature. The disappearance of structure 
with rise in temperature is also shown in the 
Original paper where the atomic density is 
plotted as a function of the distance from the 
central atom. 


‘J. A. Campbell and J. H. Hildebrand, J. Chem. Phys. 
11, 330 (1943). 
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Physics and Philosophy 


By Sir JAMEs JEANS. The Macmillan Company, New 
York, 1943. Price $2.75. 


Sir James Jeans is well known to physicists for treatises 
on mechanics and on electricity and magnetism, to astro- 
physicists for work on stellar dynamics, and in recent years 
to a wider public through a series of books discussing some 
of the implications of modern developments in physics. In 
this latest book he gives an account of the objectives of 
physics and of philosophy, and of the way in which these 
subjects have interacted as they have come in contact at 
various stages in their historical development. It is not a 
book of brilliant epigrams but is a sober and instructive 
account of the bearing that classical physics has had, and 
that modern physics may have, on various philosophical 
attitudes. 

In general, the attitude of physicists toward philosophers 
has been one of rather contemptuous superiority. It has 
often been suggested that philosophers use a multiplicity 
of words to conceal a scarcity of ideas. Indeed, Jeans’ ac- 
count of numerous philosophical arguments tends to justify 
such an attitude. It seems almost incredible that such 
childish statements and unconvincing reasoning should 
have occupied the attention of serious thinkers and have 
been the object of widespread interest and discussion. And 
yet they were and to some extent still are. Jeans suggests 
the explanation. 

The classical philosophers were concerned with the 
formulation of “‘scientific’’ proofs for the articles of religious 
faith that dominated their times, and present-day phil- 
osophers are concerned with the equivalent development 
of a philosophy of life. To this end, such problems as that 
of free will are of major importance, and are generally not 
only entirely unsolved but are not yet at all adequately 
formulated. Not only does a philosopher not know whether 
men have freedom of will or not, but he does not have a 
clear idea of what he is talking about when he mentions it. 
Physicists have excluded such questions from their pro- 
fessional consideration but have often conducted their 
practical affairs as though freedom of the will were self- 
evident. 

On the other hand, philosophers have often discussed at 
great length the discoveries of physics with only a very 
superficial knowledge of them. As Jeans points out, the 
significant part of a physical theory is its mathematical 
formulation, while as a rule the philosophers have been 
concerned with, and often only acquainted with, the verbal 
description that grows up around the mathematical core. 

After a description of some of the results of the theory of 
relativity and of various formulations of quantum me- 
chanics, Jeans proposes a rough outline of what might be 
regarded as the nature of the “real world.” In addition to 
the usual dualism of mind and matter he proposes an 
additional dualism of phenomena of which the mind can 
be conscious and a substratum entirely beyond the reach 
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of consciousness. This substratum, according to Jeans, is 
made necessary by the facts of quantum mechanics. 

One can raise a great many objections to such a formula- 
tion; but whether it or any modification of it will prove 
satisfactory must depend on whether it is of any real 
assistance in formulating philosophical problems more satis- 
factorily than has been done in the past. Certainly, the 
philosophical problems cannot be ignored. They represent 
an important part of man’s effort to adjust himself to his 
surroundings and they always demand solution. 


W. V. Houston 


Electrical Counting 


By W. B. Lewis, Pp. 144, Figs. 71, 1522 cm. Uni- 
versity Press, Cambridge, and Macmillan Company, 
New York, 1942. Price $2.50. 


This little book should be read by every serious student 
of experimental physics and should be studied by every 
nuclear physicist and electronic circuit man. No matter 
how experienced an individual may be in the field covered 
by this book he is practically certain to pick up at least one 
worthwhile idea which will be new to him. Its content is 
broader than implied by its modest title Electrical Counting 
with Special Reference to Counting Alpha and Beta Particles. 
The treatment is clear, succinct, and logical. Substantially 
every circuit problem involved in the design, construction, 
operation, and servicing of various types of counters is 
discussed. The apparatus required for such counting experi- 
ments includes common circuit elements such as voltage 
stabilizer circuits, feedback amplifiers, cathode followers, 
mixing circuits, triggered circuits, and discriminators. In- 
clusion of this material makes about three-fourths of the 
book useful and practical to all circuit designers. In addi- 
tion to this ground work, professional nuclear physicists 
will find helpful hints on the design of ionization chambers, 
the operation of discharge counters, recording oscillographs 
and message registers, the statistics of counting random dis- 
tributions, and a brief but valuable summary on alpha-ray 
energy determinations by means of range measurements. 

Dr. Lewis devotes a section to each of the fundamental 
components of detection apparatus. Each such section con- 
siders one or two basic circuit elements which experience 
has shown come closest to performing the desired functions. 
Occasionally Dr. Lewis commends a particular circuit as 
almost perfect in performance, but in every case the treat- 
ment is highly critical; the disadvantages, advantages, and 
alternatives receive proper consideration. 

Almost nothing in the book is second hand. Dr. Lewis 
‘has battled the circuits he describes, and the inventors of 
particular circuits will usually find their offspring described, 
criticized, and commended from the viewpoint of Dr. 
Lewis’ own experience with the circuits. 

Specific circuit constants are generally based on vacuum 
tubes available in England before the war. American 
workers should find that Dr. Lewis’ clear discussion of 
fundamentals will permit them better to understand their 
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own pentode-studded circuits, although the triode pre- 
dominates in British usage. No circuit developments or 
improvements made since the war began are included, but 
as a summary of principles and practice up to about 1929 
the book has undoubted permanent value. 

RosLey D. Evans 





Here and There 


Hospital X-Ray Uses Turntables 








An x-ray machine incorporating a miniature turntable 
that enables physicians to change the position of a patient 
at will during delicate operations to remove foreign bodies 
from the lungs has been developed by the Westinghouse 
Electric and Manufacturing Company for the University 
of Pennsylvania Hospital. The x-ray machine, called a 
biplane fluoroscope because it enables examinations to be 
made in two planes 





horizontal and vertical—has been 
installed in the hospital’s department of radiology. It will 
be used to locate such foreign bodies in the lungs as coins 
and safety pins. The turntable makes it possible to position 
a patient so that x-ray pictures can be made from any 
angle, even vertical, without lifting or turning him bodily. 
The result is greater safety and comfort for the patient, 
greater usefulness for the biplane fluoroscope. 


* 


Dr. Slepian gets Lamme Medal 


The Benjamin Garver Lamme Medal, awarded annually 
by the American Institute of Electrical Engineers for meri- 
torious achievements in development of electrical apparatus 
or machinery, was presented recently to Dr. Joseph Slepian, 
associate director of the Westinghouse Research Labora- 
tories. The gold medal, established in 1924 by the will of 
Mr. Lamme, chief engineer of the Westinghouse Electric 
and Manufacturing Company and an associate of Dr. 
Slepian, was presented to the winner by Harold Osborne, 
A.I.E.E. president. The ceremony took place during the 
five-day national technical meeting of the A.I.E.E. at the 
Hotel Cleveland. 


* 


Floating Lights Aid Seaplane Landings 


Because water absorbs rather than reflects light, it can- 
not be flood-lighted for airplane landings like a land run- 
way. But floating rubber ‘‘doughnuts,” carrying intense, 
battery-operated fluorescent lamps, now clearly define 
water runway areas, permitting safe landings even on the 
darkest night. For these floating beacons, Westinghouse 
engineers developed a new hot-cathode type of fluorescent 
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iamp which reduces the necessary electrical starting ‘‘pres- 
sure’ to only 45 volts—one-sixtieth that required for the 
special type of cold-cathode fluorescent lamp formerly used. 


* 


Three Harvard professors, who have recently visited 
Mexico, have been elected correspondent academicians of 


the Mexican National Academy of Sciences at a meeting 
of the Academy. Those honored were, Dr. Harlow Shapley, 
Director of Harvard College Observatory, Dr. L. C. Graton, 
Harvard geologist, and Dr. P. W. Bridgman, Harvard 
professor of physics. This is considered an important step 
in furthering closer cultural relations between Mexico and 
the United States. 





Recent Applications of Physics 








Tests for Transformers Tests made at the Sharon 
transformer works of the 
Westinghouse Electric and Manufacturing Company pro- 
vide the means for a study and reduction of transformer 
noise. Two trucks, each carrying a transformer, are wheeled 
into position in a room specially constructed to keep ex- 
terior sounds out. The doors are then closed and by means 
of instruments, located on a panel outside the vault, the 
hum emitted by each unit is measured. Truck contacts 
engage jack contacts on the inside wall of the vault; a small 
microphone picks up the sound, the relative value of which 
is read by suitable instruments. 


Increase in 
Electrical Power 


A new turbine generator of 
50,000-kilowatt capacity, suf- 
ficient to supply all the elec- 
trical needs of a city of 50,000 population, has just been 
placed in service by Consolidated Edison. Some of the parts 
for this turbine were produced by the power company in 
order to complete the machine when pressure of other war 
work halted its construction in a Westinghouse plant. Be- 
cause of the high quality of the machining operations per- 
formed by Consolidated Edison, the Westinghouse Com- 
pany subsequently let sub-contracts to its customer for the 
machining of turbine spindles, gear cases and bearings, and 
a volume of detail turbine parts for the Naval shipbuilding 
program. According to J. H. Spraggon of Westinghouse, 
this industrial teamwork provides more electrical power for 
war industries in the New York area. 


Generator Brushes 
for Aircraft 


A new method of treating 
the carbon brushes for air- 
; craft generators has been de- 
veloped by Dr. Howard M. Elsey, consulting chemist of 
the Westinghouse Research Laboratories. This develop- 
ment offsets the effects of thin air and low humidity at high 
altitudes and life of the brushes has been increased from an 
average of two hours to well over 100 hours for high 
altitude flying over 30,000 feet. These new brushes are now 
being applied in bombers and fighting planes. They were 
devised after extensive laboratory tests in which Dr. Elsey 
observed brush characteristics in a vacuum jar and in 


VOLUME 14, SEPTEMBER, 1943 


specially designed chambers simulating high altitude con- 
ditions. It is now expected that further developments, still 
on test, will provide even longer life than has already been 
made possible. 


Aircraft Stamping Dies To meet the demand for 
faster production of the mul- 
titude of blank and pierce dies needed by the aircraft and 
other industries, a 300 percent faster method of manufac- 
turing such dies has just been announced by Algoma Prod- 
ucts, 3080 E. Outer Drive, Detroit, Michigan. In addition 
the process reduces the cost 65 percent. In the Algoma 
method, dies are produced from templates of the parts 
themselves, eliminating the need for special die design. The 
construction employed simplifies alterations due to en- 
gineering changes. The usual complex, time-consuming, and 
costly milling operations are replaced by a method in which 
dies are contour-sawed from chrome molybdenum steel, 
filed, sheared to size, hardened, and welded to inexpensive 
back-up plates. One-third of the critical alloy steel re- 
quired in conventional dies is thus conserved. Individually 
replaceable pierce punches of hard drawn tool steel are 
pressed into the blank die and matching holes drilled in the 
blank punch for piercing holes in the part. 

These dies are self-stripping and self-ejecting. All punches 
and dies are treated to withstand continuous production 
service. The punch and die are kept in alignment by two 
guide pins of different diameters in one portion of the die 
and matching hole in the mating part; the different diam- 
eters prevent incorrect assembly of the dies in the press. 


Electronic Diffraction 
Camera 


Dr. Earl A. Gulbransen of 
the Westinghouse Research 
Laboratories is investigating 
the atomic structure of coatings that form on steel, alu- 
minum, and copper when these metals are exposed to air 
or corrosive chemicals. Electrons are shot through the 
three-foot long vacuum tube and bounced off a highly 
polished button of aluminum or steel on which an oxide 
coating is being built. The electrons ricochet off the faces 
of the block-like molecules that form the coating and con- 
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tinue downward at an angle to strike a strip of photo- 
graphic film. The electrons trace a pattern of black and 
white semi-circles that appear when the strip is removed 
from the camera and developed. This design is formed by 
the electrons bouncing off the different faces of the mole- 
cules in the coating. 

The electronic bullets are forced out of a small aluminum 
rod at the top of the vacuum tube by 50,000 volts. When 
the voltage is applied to the rod, a pinpoint area on the end 
of the rod gives off free electrons that speed through the 
tube at the rate of four and a half million miles an hour. 
The polished metal buttons that Dr. Gulbransen uses as 
targets are inserted six at a time in a magazine that fits into 
the side of the vacuum tube. By revolving the magazine 
he can bring any of the buttons into the field of electronic 
fire. Electric resistance wires in the magazine heat the 
buttons to accelerate the formation of coatings. Growth 
of the oxide coating can be also speeded by admitting 
measured amounts of oxygen and hydrogen into the tube. 

Dr. Gulbransen’s work with an electronic diffraction 
camera is pointing the way toward longer-wearing bearings 
and cvlinders for airplane and automobile engines, better 
tin plate in which the tin will cling more tightly to the steel 
underneath, and cheaper methods of making stainless steel. 


Walkie-Talkies The War Emergency Radio 
Service Fact Sheet gives in- 
formation regarding the small portable sending and re- 
ceiving sets carried by the defense forces. By means of 
these sets, called Walkie-Talkies, an incident officer, or a 
fire chief, or chief of any emergency unit can direct his 
squads easily and quickly from a vantage point at the 
scene of disaster. 

The new emergency radio system is of direct use to many 
strategic centers and installations in the community. Its 
importance appears in the two following representative 
cases: 

1. Industrial plants.—Calls to an industrial plant give 
air raid warning, advise of latest developments during a 
raid, confirm calls for emergency units and indicate the 
help coming, and advise in the operation of the plant’s own 
emergency forces. Calls from an industrial plant summon 
emergency medical teams, rescue units, fire and _ police 
forces, demolition squads, and if needed, decontamination 
units. Large industrial plants can use walkie-talkies to 
reach plant protection volunteers in outlying sections of 
the plant or in separate buildings or defense posts. 

2. Hospitals.—It is crucial to know during a raid exactly 
what beds are available and what operating rooms are free 
in the casualty receiving hospitals of a community, and to 
’ direct casualties rapidly to available facilities. Ambulances 
can be loaded and dispatched effectively because the con- 
trol center is in communication both with hospitals and 
with the incident officers and incident medical officers. 
Mobile medical teams can be directed from point to point 
without returning to their bases. If the hospital facilities 
of an area become overtaxed, the radio channel to the 
district headquarters can arrange for reinforcing hospital 
facilities. 
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WALKIE-TALKIE GOES TO SEA 


Although developed to ride a soldier’s back in steaming 
jungles, desert wastes, or arctic drifts, the “‘walkie-talkie,”’ 
one-man radio station built by Westinghouse for the Army, 
now has gone to sea. Its ten-mile range is short enough so 
that ships in convoy can transmit voice messages from 
vessel to vessel without using the “‘big’’ radio, which might 
give away ship location to enemy listeners. 





Innovations in Instruments 








New and Rare Instruments 


Information concerning the following offers and requests 
for rare instruments that can be sold, loaned, or leased for 
essential war research work can be obtained from D. H. 
Killeffer, Chairman of the Committee on Location of New 
and Rare Instruments, 60 East 42nd Street, New York 17, 
New York. 


Instruments Offered: 


L & N Thermionic Amplifier; Ammeters—Whitney and 
Jewel; Various Balances; Schmidt & Haensch Colorimeter; 
Platinum Calorimeter; Weston D. C. Electric Meters— 
Voltmeters, Ammeters, Galvanometers; Two-Circle Goni- 
ometer (two loan for war use); Microscopes—Bausch & 
Lomb Centrifuge, Zeiss Binocular, Leitz (suitable for 
photomicrography); Poeller Optical Wedge Pyrometer; 
4L&N Portable Potentiometers; L & N Type K Potenti- 
ometers; Polarizer and Analyzer for microscope; Abbé 
Refractometer; R. Fuess Reflectometer; Spectrographs, 
Spectroscopes, Spectrophotometers; Saccharimeters; Say- 
bolt Universal Viscosimeter. 





Instruments Requested: 


Western Electric Electrometer; Flexaform Set; Amsler 
No. 4 Integrator; Micromanipulators. 


Student’s Kelvin Bridge 


A new illustrated 8-page booklet issued by the Leeds & 
Northrup Company describes the Student’s Kelvin Bridge 
[Catalog E-53(1)], an instrument which provides the fea- 
tures which teachers have found most valuable in imparting 
an understanding of the Kelvin bridge method. Although 
intended primarily for instruction, the bridge is said to be 
suitable for general use in the measurement of resistance 
from 0.1 to 0.001 ohm, when only moderate accuracy is 
required. Its design is simple, and affords a clear illustration 
of the fundamental principles of the method. Provision is 
made for inserting extra resistances in either of the ratio 
circuits or in the yoke, so that the effects of resistance in 
these parts of the bridge can be demonstrated. Copy of this 
booklet sent upon request by Leeds & Northrup Company, 
4934 Stenton Avenue, Philadelphia, Pennsylvania. 
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The Approximate Mathematical Methods of Applied Physics as Exemplified 
by Application to Saint-Venant’s Torsion Problem 


THOMAS JAMES HIGGINS 
Illinois Institute of Technology, Chicago, Illinois 


(Received July 7, 1943) 


HIS paper is one of a series of three papers 
devoted to exposition of the known solu- 
tions of Saint-Venant’s torsion problem for a 
prism of uniform cross section. The first,! pub- 
lished some time previously, comprises account 
of the historical development of the fundamental 
mathematical and physical theory of the torsion 
problem and of the solutions that have been 
effected by exact analytic methods. The third,? 
to be published shortly, covers the solutions that 
have been effected by empiric and by experi- 
mental methods. The present paper, the second 
of the series, embraces those solutions that have 
been effected by various approximate methods. 
The more important of these methods com- 
prise three categories: variational, arithmetic, 
and graphic—according as they stem from the 
theory of the calculus of variations, from the 
theory of linear difference equations, or from the 
use of graphic analysis. In that they afford, in 
divers branches of elastic theory, solutions of 
problems mathematically intractable to exact 
analysis (which category includes, for example, 
most of the important practical torsion problems 
associated with the stress analysis of high speed 
airplane, engine, and tool designs) these methods 
are of considerable technical importance. Yet, 
despite this significance, no considerable exposi- 
tion of these approximate methods is to be found 
in standard reference texts on elastic theory or 
applied mathematics; nor for that matter, is other 
than limited mention* made therein of sources, 
either original or secondary, that contain such 
exposition. 
A host of such sources, however, is mentioned in 


*In fact, none whatsoever of most of the methods dis- 
cussed in this paper. 
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this paper. And in that these furnish a connected 
and logical introduction to a number of the more 
powerful approximate mathematical methods 
available for effecting actual numerical solutions 
of two-dimensional equilibrium problems in all 
branches of applied physics, especial interest may 
attend the following inclusive, and fully docu- 
mented, account of the application of these 
approximate methods to Saint-Venant’s torsion 
problem.** 


I. VARIATIONAL METHODS?+ 
1. Ritz’s Method and Modifications 


Thirty-five years ago (1908) the Swiss physi- 
cist Walther Ritz published a memoir* wherein 
he propounds a new method of solving problems 
of mathematical physics; and in this and sub- 
sequent memoirs,*® he effects by this method 
approximate solutions of problems which had 
long defied the efforts of analysts. Among the 
classical problems treated in the first memoir is 
the Dirichlet problem for a singly-connected 
region. Inasmuch as the torsion problem for a 


** For a qualitative account of the torsion problem see 
reference 1; for accounts of the general mathematical 
theory see: S. Timoshenko, Theory of Elasticity (McGraw- 
Hill Book Company, New York, 1934); I. S. Sokolnikoff, 
Mathematical Theory of Elasticity (Brown University, 
Providence, 1941); A. E. H. Love, The Mathematical Theory 
of Elasticity (University Press, Cambridge, 1934), fourth 
edition. 

+ The notation and terminology of this paper are those 
of Love: » denotes the modulus of rigidity; r the angle of 
twist per unit length; and the other symbols are defined 
as introduced. To avoid redundance symbols are used in 
the text in lieu of the corresponding phrase. Again, inas- 
much as determination of either the torsional moment M 
or the torsional rigidity C is tantamount to determination 
of the other, for (M=7rC), mention is made only of M— 
though in some of the papers discussed it is the expression 
for C, and not M, that is given. 
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“‘singly-connected” prismff is interpretable as a 
first Dirichlet problem (namely, as the search for 
the conjugate torsion function y that satisfies 
V*y=0 over the right cross section S of the 
twisted prism and satisfies y = (x?+-y?)/2+<a con- 
stant on the boundary s of S) the solution for a 
particular cross section is afforded by Ritz’s 
method—sometimes termed the Ritz-Rayleigh 
method, Lord Rayleigh*® having pointed out that 
Ritz’s procedure has much in common with 
approximate methods used earlier by him. 
Application of Ritz’s method to a problem in 
elasticity hinges on expressing the problem as one 
in the calculus of variations, a transformation 
usually effected with ease inasmuch as the de- 


finitive partial differential equations of most. 


physical problems can be formulated from one 
of the well-known ‘“‘principles’’ of physics—for 
example, Hamilton's, or Maupertuis’, or Cas- 
tigliano’s principle. Transformed, the problem 
reduces to determination of a function which 


ae 


minimizes a certain integral expression subject 
to specified boundary conditions. The torsion 
problem, for instance, can be interpreted as the 
search for the stress function ¥[ =y— (x?+ y*)/2] 
that minimizes the integral expression 


UW) =(ur?/2) f grad wdS 


for the elastic energy per unit length UL=7rM/2] 
subject to the conditions: (i) the twisting couple 


M=2ur { vas 
S 


is a prescribed constant; (ii) ¥=a constant (com- 
monly taken as zero) on s. 

An alternative, somewhat more convenient 
statement of the problem is to find WV that 
minimizes the integral expression* 


tt The torsion problem for a prism of multiply-connected 
cross section is also solvable by Ritz’s method (see reference 
15); however, the boundary conditions—more complicated 
than for the “‘singly-connected” prism—are such that 
characterization as a first “Dirichlet problem” is hardly 
proper. 

* Note that by virtue of the quantities defined in the 
foregoing paragraph J(¥) = U(¥)— Mr; and since M is a 
prescribed constant, J=—1rM/2=-—U, wherein I and U 
designate the minimums of J(W) and U(¥). 
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I(¥) = (ur?/2) f (grad v—4)dS 
s 


subject to the boundary condition ¥=0 on s. 
And yet other forms of the integral expression to 
be minimized by the same W are to be found. 

The tie between the two methods of inter- 
preting the torsion problem—as a variational 
problem or as a boundary-value problem in 
partial differential equations—lies in this: The 
W which solves the variational problem neces- 
sarily satisfies** a characteristic partial differ- 
ential equation (also termed the Euler equation) ; 
this equation is precisely the partial differential 
equation occurring in the boundary-value prob- 
lem; and inasmuch as the boundary condition is 
the same in both problems—necessarily so,*** 
in fact—it follows that this V is also that which 
solves the boundary-value problem. 

In general, Ritz’s method does not yield the 
function W that exactly minimizes the integral 
expression J(W): Rather, it vields an mth ap- 
proximation ¥, to ¥Y—where Y, is expressed as 
a linear combination of » functions of a certain 
type individually satisfying the boundary condi- 
tions—and an mth approximation J,(¥,) to J. 
But often the approximation to the true minimum 
I can be effected to any desired degree in a 
comparatively simple manner; the approxima- 
tion to the true minimum constantly improves 
as the terms of the assumed series increase in 
number; and if an appropriate selection of func- 
tions has been made—commonly, this selection 
is determined by experience, by the particular 
coordinate system employed, and by the fact 
that the integrand stemming from substitution of 
VY, must be mathematically tractable—the ap- 
proximation is often quite good when the number 
of terms in W, is small. In general, the security 
for ¥, and J, being actual approximations to ¥ 
and IJ is ensured by satisfaction of the conditions 
that stem from setting the problem as one in the 
calculus of variations. Specifically, exact condi- 
tions on W sufficient to ensure the approach of 
Vv, and J, to VW and J have been the subject of 








** See A. R. Forsyth, Calculus of Variations (University 
Press, Cambridge, 1927), p. 551. 

*** See reference 15, pp. 4-7 and the extension of this 
discussion in reference 35, Chapter 3. Also, H. Bateman, 
Partial Differential Equations of Mathematical Physics 
(University Press, Cambridge, 1932), p. 157. 
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research by Ritz, by Trefftz, and by others.f 
Provided these sufficient conditions are satisfied, 
they have shown that 


lim V,=V; 


lim J,(¥,)=Z; and I[,(¥,)>T. 
nD nD 
Now, though Ritz’s method yields a solution 
where classical methods fail, and though by 
virtue of the inequality J, >J this solution is an 
upper limit to the exact solution, it commonly 
ensues that the degree of approximation to the 
exact solution is unknown. In appreciation of 
this difficulty Trefftz has developed a method, 
akin to Ritz’s, but so designed that J,’(=In) 
approaches the true minimum J from below. As 
with Ritz’s method, the approximations to V 
and to J can be carried to any degree of ac- 
curacy ; and providing W,, satisfies certain known 
conditions, 


lim ¥,=V; 


nD 


lim I,’(V¥,) =I; 


nD 


but now I,’(¥,) <J. Thus, by application of both 
Ritz’s and Trefftz’s methods the true minimum 
is bounded in known: fashion, J,’<J<J,; and 
the sign and magnitude of the maximum possible 
deviation of J,, J,’, or their arithmetic mean 
(possibly a closer approximation) from J is 
determinate. 

Trefftz has developed the theory of his method 
in two memoirs. One,’ essentially mathematical, 
is devoted to elaboration of the general theory of 
his method. The other,* essentially expository, 
comprises solution of the torsion problem for 
certain cross sections: ¥, and M, of a [prism of ] 
square [cross section ] are obtained by both his 
method and Ritz’s and the two solutions com- 
pared in various fashions; a device for extending 
both methods to provide solutions for cross 
sections having re-entrant angles is outlined ; by 
both methods M, is effected for an angle-iron 
with identical, finite flanges. 

Consideration of solutions effected by both 
methods indicates that, in general, certain ad- 
vantages attend the use of Trefftz’s rather than 


} For a number of references to these purely formal in- 
vestigations see the bibliography of reference 15; see also 
N. Kryloff and J. Tamarkin, Bull. Acad. Sci. Russ. [6] 12, 
69-88 (1918); E. Trefftz, Handbuch der Physik (J. Springer, 
Berlin, 1928), Vol. 6, pp. 130-132. 
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Ritz’s method. Thus—as Trefftz professes 
although the approximating functions used in 
the latter need not necessarily satisfy the partial 
differential equation defining the problem, they 
must satisfy the boundary conditions—not 
always a simple thing to ensure. Contrariwise, 
Trefftz’s method—to be characterized as a 
generalization of the classical method of develop- 
ment in terms of particular solutions—requires 
only a choice of functions that satisfy the partial 
differential equation, commonly a simpler de- 
mand. Secondly, the choice of functions in 
Trefftz’s method is usually such that determina- 
tion of the maximum error in the solution is 
easier. Finally, proof of convergence of the 
approximating function to the true minimum is 
generally simpler. One considerable advantage, 
though, attends Ritz’s method: ordinarily, the 
approximation converges somewhat faster than 
in Trefftz’s method. Thus, Holl and Anderson,?® 
using both methods to obtain values of M, for 
the cross section bounded by arcs of two confocal 
parabolas,* find that though both values are in 
good agreement with the known exact solution 
the value afforded by Ritz’s method is in closer 
agreement. 

A second method of obtaining a lower bound 
to the minimum function J solving a given 
problem, one in some ways more general than 
that of Trefftz, has been outlined by Friedrichs."® 
It consists of associating with the given minimum 
problem the auxiliary problem of finding the 
maximum of a second integral, obtained from 
the first by a simple transformation and of such 
nature that its maximum is identical with the 
minimum of the first integral. As, obviously, the 
negative of the maximum sought is a minimum, 
it can be obtained in a fashion akin to that of 
Ritz. The actual procedure can differ in that if 
only the approximation J, to the extremum J is 
sought, and not, as well, the approximation ¥, 
to VW producing J, the approximating functions 
need not necessarily satisfy the boundary condi- 
tions as they must in the Ritz process. On the 
other hand, convergence to the minimal solution 
is somewhat slower. 

The explicit theory essential to application of 


* Statement of an mth approximation implies determina- 
tion of the preceding »—1 approximations. 
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Friedrichs’ method has been developed by 
Friedrichs himself and by Basu." The latter 
treats in detail an important topic barely touched 
by Friedrichs: namely, determination of the 
[torsion] function ¢, the complex conjugate 
of ¥. If WV, is obtained by Ritz’s method, y, 
follows immediately from the relationship 
y= V+(x°+y")/2. But though ¢ is the complex 
conjugate of y, ¢, cannot be obtained from y, 
by the usual methods of function theory: for 
the terminated series expressing V, does not, in 
general, satisfy the differential equation of the 
problem V?¥= —2, even though the conditions 
on ¥, are such that the function represented by 
the limit of the series as n—>~ does satisfy the 
equation term by term (an interpretation of 
this latter mathematical phenomenon has been 
advanced by Goldsbrough") ; whence it follows 
that y, obtained from ¥,, is not a harmonic func- 
tion. Friedrichs merely indicates how deter- 
mination of @ conjugate to y can be expressed 
as a variational problem. Basu develops in detail 
the required analysis to the end that ¢ is ob- 
tained through determination by Ritz’s method 
of the function J maximizing a certain double 
integral expression. Additionally, the approxi- 
mation to this integral expression furnishes the 
desired lower bound on J, or J, <J<TJ,. Further, 
by virtue of the identity J= — 7M/2, bounds are 
provided on M: that is, J,<7.M/2<J;,. 

The details of the analysis contained in Basu’s 
memoirs are exemplified by application to the 
torsion problem of the square. In the first 
memoir'! WV» is effected by Ritz’s method; M2 by 
both Ritz’s and Friedrichs’ methods. In the 
second memoir" ¢2 is obtained by Ritz’s method ; 
and from @2 another value of M2. Close agree- 
ment of this last value of M2 with the known 
value is taken as evidence that the corresponding 
expression for ¢2 is sufficiently accurate for most 
purposes; i.e., to give the longitudinal displace- 
ments of points of the initially plane cross 
section. 

At this point we must mention a certain defect 
in the methods discussed above. Application of 
Ritz’s and of either Trefftz’s or Friedrichs’ 
method enables estimation of the error in an 
approximation to the minimum I of the integral 
expression J(¥); for these methods furnish 
upper and lower bounds on J. Often, however, 
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what is desired, yet usually cannot be obtained 
inasmuch as the mathematical technique is 
lacking, is knowledge of the error in the approxi- 
mations afforded by WV, to V or to its derivatives. 
In the torsion problem, for instance, the longi- 
tudinal displacement and the stress distribution, 
respectively furnished by ¢ and by the deriva- 
tives of WV, are of greater immediate value in 
practical design than is either U or M afforded 
by J. But though, as just remarked, a general 
technique is not, available for calculating the 
error in V,, it happens that a specialized tech- 
nique applicable to the torsion problem is known. 

In a paper devoted to exposition of his 
“thickness parameter’? method, explained in 
Section II below, the English aeronautical en- 
gineer W. J. Duncan" outlines a method* of 
calculating the maximum deviation between the 
values of a function WV, defined by the conditions 
that it satisfy V?¥ = —2 over S and vanish on s, 
and those of an approximation V,. Thus, if é’ 
and ¢” are the maximum and minimum values 
within S of a function ¢ defined by V?¥,+2=e, 
then (1—e”/2)=W,/W=(1—€'/2) ; and conjunc- 
tion of this knowledge with the maximum values 
of € on the boundary enables determination of the 
maximum deviation of ¥, from W at any point 
within or on the boundary of S. 

This procedure, however, does not furnish—at 
least directly—information on the derivatives of 
V,; nor on ¢, the conjugate of 


yL=¥+(x?+ 9") /2]. 


Lacking means to obtain such information, it is 
tempting to conjecture that if J is bounded 
closely by J, and J,, the corresponding expres- 
sions for ¢,, V,, and the derivatives of the latter 
yield values in good agreement with the exact 
values. Though such inference is as yet unsub- 
stantiated by rigorous analysis, empirical support 
might derive from comparison of values so ob- 
tained with exact values calculated from known 
rigorous solutions (see reference 1 for discussion 
of more than a 100 such solutions). 

In a recent paper whereof several pages are 
devoted to a qualitative discussion of the funda- 
mental concepts of the Ritz process and of its 


* Obviously, this method is directly extensible to the 
more general equation, V2¥ =an arbitrary constant. 
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merits and defects as regards actual use 
discussion 


which 
complements foregoing remarks— 
Courant,'® illustrative of his remarks on the 
solution of the torsion problem for multiply- 
connected cross sections, obtains V. and M, for 
a hollow square. This, to the writer’s knowledge, 
is the only published treatment of the numerical 
solution of a multiply-connected cross section by 
Ritz’s method. 

Finally, we note two other applications of the 
Ritz analysis to the torsion problem. Approxi- 
mating the normal cross section of a propeller 
blade by the area bounded by arcs of two non- 
orthogonal parabolas, Leibenzon'® obtains 
and a first approximation 7; to the maximum 
shearing stress T,,; for an area bounded by one 
of the parabolas and by a line segment normal 
to its axis, YW. and 7,2; and for a streamline 
section akin to an airfoil, W2, Ms, and T,,2 (how- 
ever, Tm» is obtained from ¥,, by differentiation, 
in which connection see the second preceding 
paragraph). 

Duncan, Ellis, and Scruton!’ utilize Ritz’s 
method to investigate the torsion problem for the 
general isosceles triangle, the exact analytic solu- 
tion of which is as yet unknown, though certain 
special cases have been so solved. V, and M, are 
effected. It is found of the former that it is exact 
for the special case of the equilateral triangle and 
is a good approximation for isosceles triangles with 
small included angles, except, possibly, near the 
base of the triangle; and of the latter that it is 
“almost certainly amply accurate for practical 
purposes.”’ 


2. Galerkin’s Method 


In 1915 Galerkin'® advanced a variational 
process which has found considerable use in 
divers fields of applied mechanics, especially in 
aircraft stress analysis. Galerkin’s method has 
much in common with Ritz’s process: indeed, 
Biezeno,'!* demonstrating that through analytic 
transformation the general equations of one can be 
derived from those of the other, maintains that 
in essence they are identical. Elsewhere, dis- 
cussing the relative merits of Ritz’s, of Galerkin’s, 
and of a third method advanced by Koch?® and 
himself, Biezeno*! indicates a process that en- 
compasses the just-named as special cases. More 
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inclusively, Gross? has shown that these pro- 
cedures and others are derivable as special cases 
of a very general functional method. 

In Galerkin’s method an mth approximation to 
the rigorous solution of a one-dimensional 
problem is taken in the form 


4 cifi(x), 


wherein the c; are arbitrary constants and the f; 
a selected set of linearly independent functions 
individually satisfying the boundary conditions. 
Substitution of this expression in the defining 
differential equation (arranged with right-hand 
member zero) yields ‘‘a remainder”’ or ‘‘error”’ 
e, a linear function of the ” arbitrary constants. 
Then multiplication of « by each of a set of n 
arbitrarily chosen, linearly independent func- 
tions termed multipliers; integration of each of 
the products thus formed over the region S of 
the problem; equation to zero of the resulting 
expressions; and solution of the resulting set of 
equations yield a preferred set of values for the 
arbitrary constants. 

From the standpoint of mechanics e is inter- 
pretable as a generalized force, the arbitrary 
constants as a set of generalized coordinates, the 
appropriately chosen multipliers as virtual dis- 
placements corresponding to increments of each 
of these generalized coordinates, the integrated 
products as weighted means of ¢, and the vanish- 
ing of these weighted means as the vanishing of 
the virtual work due to the corresponding dis- 
placements. In light of this mechanistic inter- 
pretation it follows that the degree of accuracy 
obtained can be increased indefinitely by in- 
creasing the number of functions employed ; that, 
in general, exact solution of an m-dimensional 
problem requires an m-fold infinite set of func- 
tions fi(x1, «++, Xm); but that if the functions and 
the multipliers are chosen appropriately, it can 
be anticipated that a good approximation stems 
from use of only a few functions. 

Now in terms of the familiar membrane anal- 
ogy to the torsion problem ¥ is proportional to 
the vertical displacement of the membrane and 
€ is proportional to the applied pressure not 
balanced by the tensions in the membrane. Ac- 
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cordingly, we represent the stress function by 


n 
V=)>  cifi(x, y), 
os 
these functions f; appropriately chosen and satis- 
fying the boundary condition ¥=0 on s; sub- 
stitute V in V?¥+2=0; obtain the “‘error,” 


n 
e=)> Dici, 
i=1 

wherein the ); are known constants; choose as 
multipliers the f;; integrate each of the products 
ef (x, y) over S; equate to zero each of the result- 
ing expressions—hence, impose the condition 
that the work done by the unbalanced load in a 
virtual displacement of the membrane propor- 
tional to f(x, y) be zero; and solve the resulting 
set of algebraic equations for the values c¢;. 

The most inclusive treatment* in English of 
Galerkin’s method is encompassed in the writings 
of Duncan.*-*7 Illustratively, Duncan*:?5 applies 
this theory to a variety of problems, among them 
the torsion problem. Thus, for an airfoil bounded 
by the loop of a bisymmetrical cubic oval ¥, 
and Ms are effected and compared with cor- 
responding values obtained earlier by Duncan 
by his “thickness parameter’’ method. Analysis 
of a bisymmetrical parabolic lenticular cross 
section is carried out by both methods: VW and MW 
are calculated, respectively, to the 4th and 6th 
powers of the thickness parameter ; the latter is 
found to be in close agreement with MM, effected 
by Galerkin’s method. 3; is determined for a 
rectangle; values of MW; for squares are in good 
agreement with values obtained from Saint- 
Venant’s exact solution. Finally, ¥; and M, are 
obtained for an arbitrary area possessing at least 
one axis of symmetry ; corroboratively, M, yields 
the known exact expression for the ellipse. 


3. Other Variational Methods 


A variational method proposed by Poritsky?$ 
is of considerable usefulness. As exemplified by 
the torsion problem ¥ is taken to be of the form 


n 
WV = z=. Xj Y;, 
i=1 
* See also a very brief discussion in E. Keller, Mathe- 
matics of Modern Engineering (John Wiley and Sons, New 
York, 1942), Vol. 2, pp. 288-292. 
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wherein the Y; are assumed functions of y alone 
and are such that W satisfies the boundary con- 
dition ¥=0 on s, and the X; are unknown func- 
tions of x alone. Then substitution of V in I(¥) 
and imposition of the variational condition that 
6I1(v) = Y,6X; vanishes for all 6X; yields a set of 
n differential equations, each containing but one 
X;. Whence the set of X; is determined by 
finding a set of solutions that satisfy boundary 
conditions determined in part by the assumed 
form of the Y; and in part by the original bound- 
ary condition. 

In this fashion Poritsky obtains ¥,, and from it 
M, and T,,;, for a certain rectangle. These values 
are in excellent agreement with corresponding 
values calculated from Saint-Venant’s exact 
solution. In similar fashion WV, is obtained for the 
ellipse and for the rhomboid. 

In solution of the torsion problem, or yet other 
two-dimensional problems pertinent to a given 
area, certain modifications of the above pro- 
cedure may prove advantageous: say, inter- 
change of the values of X; and Y;; or choice of 
Y; as a function of both x and y; or utilization 
of some other permissible variant. 

Yet another variational method is described 
by Hovgaard.?® By conjunction of Castigliano’s 
principle and the use of Lagrange’s multipliers 
general variational equations are derived for a 
prism constrained to simultaneous flexure and 
torsion. Though the details of solution are not 
given, Hovgaard states that these equations have 
been used to solve the torsion problem for the 
ellipse, quasi-ellipse, equilateral(?) triangle, and 
rectangle. 


II. DUNCAN’S METHOD OF “THICKNESS 
PARAMETER” 


Of especial interest in aircraft stress analysis is 
solution of the torsion problem for uni- or bisym- 
metrical cylinders of cross section such that all 
values of the half-thickness ¢ are small compared 
with the length of the longest (major) axis of 
symmetry. Particular interest, therefore, accrues 
to methods especially adapted to the solution of 
such thin cylinders. Possibly the most effective 
is Duncan’s'* method of “thickness parameter.”’ 

Directing the x axis of a Cartesian system 
along the major axis of the cross section, the 
boundary equation is written y=t= 67, wherein 
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6 is the thickness parameter independent of x 
and T is an assigned function of x. Now, by 
virtue of symmetry about the x axis only even 
powers of ¢ and y can appear in WV; further, by 
virtue of the boundary condition on ¥ it vanishes 
on s. These considerations suggest expression of 
Vv in the form 


oc 
V=P?—y? +) (t"—y™)fr(x), t=OT, 
i=1 

wherein the functions f,(x) are to be determined 
from the condition that WV satisfies V?¥=—2. 
The f,(x) determined—discussion of the actual 
procedure is out of place here—W appears as a 
convergent series in powers of 6? with coefficients 
which are functions of x. Commonly, W is deter- 
mined only to some specified power (say, the kth) 
of 6; hence, ¥.~ WV. With Y; determined, the 
torsional moment and the shearing stress are 
obtained from the known relationships in which, 
however, WV; replaces ¥—though no proof of 
validity of this procedure is advanced. Nor is 
proof given that 





lim ¥,.=V: 


kn 
though inasmuch as the procedure (Duncan’s) 
discussed in Section [:2 affords calculation of 
the maximum error in ¥;, this defect is of aca- 
demic rather than of practical significance. 

Illustrative of his analysis Duncan obtains 
the exact form of W for the ellipse, equilateral 
triangle, thin isosceles triangle with two infinite 
sides, and curvilinear isosceles triangle comprised 
of two linear sides and base a hyperbolic arc 
subtending 45°. An area bounded by the loop 
of a cubic curve is treated at length: Ws, Myo, and 
S; are obtained; a detailed calculation of maxi- 
mum error in Wy is set out. 

In a succeeding paper Duncan*® extends the 
method of thickness parameter to doubly-con- 
nected and to asymmetrical cross sections. A 
very complete discussion of the general theory is 
illustrated by determination of approximate 
values of W and M for the slit and the complete 
annuli, and for a cross section bounded by 
concentric similar ellipses. 


Ill. FINITE DIFFERENCE METHODS 


The theory of finite differences affords power- 
ful methods of effecting numerical solutions to 


VOLUME 14, SEPTEMBER, 1943 


boundary value problems which are mathe- 
matically intractable to attack by exact analytic 
methods. In one fashion the general scheme of 
solution involves (i) replacement of the partial 
differential equation and the boundary conditions 
defining the problem by an analogous partial 
difference equation and equivalent boundary 
conditions; (ii) effectance of a sequence of ap- 
proximate numerical solutions of the difference 
problem; (iii) application of a proper test to 
determine the convergence of this sequence to 
the exact solution; (iv) determination of the 
limits of error in the final approximation. In yet 
another fashion the problem can be set up as a 
variational problem involving minimization of 
an integral expression, this latter paralleled by 
its finite difference analogue, and approximate 
solution of this last effected in an appropriate 
fashion.* 

Though solution of boundary value problem 
by application of finite difference methods is an 
idea far from new, the pertinent theory is, 
seemingly, but littke known among applied 
physicists, and even less so among engineers. In 
part, at least, this state of affairs is attributable 
to the lack of a suitable reference text in English. 
Thus, though L. M. Milne-Thomson’s The Cal- 
culus of Finite Differences covers the formal 
mathematical theory of finite differences, an 
equally authoritative text on applied finite-dif- 
ferences (aside from actuarial use) is yet to be 
written. Perforce, then, the earnest worker must 
turn to the periodical literature on applied finite 
differences. Hereof, a logical course of study 
might well encompass reading of the work of, 
say, Richardson,*"* Liebmann,** Wolf,* and 
Courant**—who present general theory from a 
semi-applied point of view—and subsequent 
study of a sequence of papers wherein develop- 
ment of theory is illustrated by application to 
one of the classical problems of the clamped 
plate, the membrane, or the torsion problem.f 

* A good illustration of this second kind of analysis is 
contained in S. J. Newing, Phil. Mag. 32, 33-49 (1941). 
Incidentally, the analysis of reference 31 is extended herein. 

7 Of this literature no considerable published bibli- 
ography is known to the writer. Of his own inclusive bibli- 
ographies on the three problems mentioned, that pertinent 
to the torsion of uniform prisms is presented in this paper; 
that to non-uniform but axially symmetrical prisms is to 
be presented in a subsequent paper. 


In a recent paper (reference 15) wherein the funda- 
mental concepts of certain finite difference procedures are 
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To the writer’s knowledge the first solution of 
the torsion problem by finite differences was 
effected by Runge,** who determined and plotted 
the stress lines in a cross-shaped area comprised 
of five squares. Extension of Runge’s analysis 
to multiply-connected areas is discussed by 
Neményi,*’ though he advances no explicit solu- 
tions of the torsion problem. 

In 1918, in a paper frequently mentioned in 
the subsequent literature, Liebmann* suggested 
certain modifications in the methods then used 
for solving physical problems by finite differences. 
Subsequently, these suggestions were reworked, 
improved, and extended by many writers, to the 
end that comparison of their papers reveals such 
overlapping that it is a matter of some difficulty 
to assign specific improvements in technique to 
given writers. 

Gerschgorin has written a series of papers, for 
the most part in Russian, wherein linear differ- 
ence equations are applied to effect the solution 
of problems involving Laplace’s and Poisson’s 
equations. Fortunately for the linguistically 
limited, one** of the more important papers is 
available in German. Therein, methods for cal- 
culating the error pertinent to the approximate 
solutions of a linear partial differential equation 
of elliptical type (the ubiquitous Laplace and 
Poisson equations are of this type) by difference 
processes are discussed. The theory developed in 
this paper is exemplified by approximation of V 
for a specifically dimensioned rectangle, deter- 
mination of the error at a selected point of this 
rectangle, and check of this latter through com- 
parison with corresponding values determined 
from the known solution. Collatz*® has extended 
the general analysis of this paper. 

Schwalbe,*® also, has studied solutions of the 
rectangle by finite difference methods. Numerical 
values for M and T,, are found to be in good 
agreement with corresponding values obtained 
from the exact solution. Additionally, finite dif- 
ference analyses are effected for the torsional 
moments of a sharp-cornered channel of uniform 





briefly but interestingly discussed, Courant mentions a 
number of papers on finite differences. Many of these 
papers, though, are on purely formal points of the theory 
and contain no illustrative examples; however, taken as a 
third unit in the above scheme of study, they can be read 
with profit by those possessing the required mathematical 
background. 
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thickness and an I beam with sloping sides. 
Numerical values calculated therefrom are com- 
pared with corresponding values obtained both 
from experiment and from various analytic 
solutions. 

Certain improvements in the original Lieb- 
mann process are professed in a series of loosely 
connected papers by Shortley and Weller ;*' 
Weller, Shortley, and Fried; Frocht and 
Leven ;* and Poritsky, Snively, and Wylie.“ 
Essentially, these papers are concerned with the 
development of methods that minimize the nu- 
merical labor incident to the solution of a given 
problem. Fuller discussion here is not feasible. 
The theory developed in the second of these 
papers is illuminated by solution of M and T 
for two squares with rounded corners—the sides 
are equal, the radii different. 

In the course of the past few years Southwell*® 
and a group of fellow workers have applied a 
numerical method, the ‘“‘relaxation method,” to 
effect the solutions of divers types of physical 
problems. Originally devised by Southwell for 
the determination of stresses in redundant frame- 
works—in which connection it has much in 
common with Hardy Cross’ **® method of distrib- 
uted moments—it was subsequently developed 
to the end that in an investigation of the general 
theory of the application of relaxation methods 
to linear systems, Temple‘? remarks that ‘‘almost 
all linear systems of equations, algebraic, integral, 
or differential, can be brought within the scope 
of relaxation methods, which seem to constitute 
one of the most powerful methods of computa- 
tion in mathematical physics and engineering.” 
Investigations of necessary and sufficient con- 
ditions ensuring that the approximations stem- 
ming from relaxation methods converge to the 
exact solutions are to be found in papers by 
Black and Southwell** and by Temple.” 

Christopherson and Southwell*® have con- 
joined the relaxation process with a finite dif- 
ference procedure indicated earlier by Thom*® to 
effect the solution of the torsion problem for an 
equilateral triangle pierced by three symmetri- 
cally located circular holes. V is approximated and 
the stress lines plotted. Preliminary to the solu- 
tion of this problem and indicative of the best 
procedures to be followed in solving it, ¥ of the 
equilateral triangle is approximated. 
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Earlier, Orr®™ had used an extension of Thom’s 
method to effect approximations to M and T, 
for a serrated annulus and for a hollow square. 
Additionally, approximations to M for a number 
of standard structural sections are determined, 
and then compared with experimental values. 

Finite difference analysis developed by Suna- 
tani and Negoro® has been employed by the 
latter® to study the stress distribution and to 
determine T,, in squares pierced by a centrally 
located circular hole. Values of the latter are in 
general agreement with corresponding . values 
determined by experiment. Additionally, an 
estimation of accuracy is afforded by comparison 
of approximate values calculated for T,, in the 
annulus with values calculated from the known 
solution ; and by comparing values of 7, for the 
pierced squares with corresponding values ob- 
tained by experiment. 


IV. GRAPHIC METHODS 


Industrial stress analysis frequently requires 
solution of the Dirichlet problem for an area 
bounded by a curve unamenable to simple mathe- 
matical expression. In this case recourse to some 
form of experimental or graphic solution may 
prove profitable. In this connection special inter- 
est attaches to Bairstow’s™ ingenious analytic- 
graphic method based on a theorem due to 
Green: “. .. any irrotational acyclic fluid 
motion can be reproduced by an appropriate 
choice of simple sources round the fixed bound- 
ary.’’ As applied to the torsion problem, a process 
of repeated contour integration of the known 
boundary values of ¥[=(x?+y*)/2 on s] yields 
the strength x of these boundary sources. 
Theoretically, rigorous determination of x re- 
quires, in general, an infinite number of repeated 
integrations; practically, a good approximation 
Xn to x is commonly afforded by calculation of 
only a few terms of the sequence. x, determined, 
the torsion functions are found by appropriate 
integrations or differentiations, effected either 
analytically or graphically. 

Illustrative of theory Bairstow, using graphic 
integration, obtains x2 for the ellipse, xo for the 
eccentric annulus, and x; for a cross section akin 
to that of a propeller blade. For none of these 
configurations does Bairstow prove the con- 
vergence of x, to x; nor is a proof of convergence 
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for the general case advanced. In the same 
paper, however, Perry (writing separately) 
through analytic evaluation of the integrals 
occurring in the solution of the ellipse and the 
eccentric annulus shows that in these cases, at 
least, the approximate solutions converge to the 
known exact solutions. 

In a subsequent paper Bairstow and Pippard®® 
advance an inclusive account of the theory of 
this method : Each step is justified rigorously and 
is interpreted physically. Illustratively, the 
scheme of solution and the pertinent mechanical 
details of the graphic analysis are delineated 
through solution of numerical values of x, ¥, ¢, 
M, and T for an ellipse. Excellent agreement of 
the calculated numerical values for the shear 
stress on the boundary with those calculated 
from the exact analytic solution are interpreted 
as justification of the accuracy of the graphic 
procedure. Less detailed treatments are given of 
several configurations of practical interest: a 
slotted circle, a slotted eccentric annulus, and a 
series of squares with rounded corners. These 
shapes approximate the cross sections of standard 
types of shafts. 

We note here that a decade before Bairstow 
advanced his method Runge*®® outlined a method 
of solving a Dirichlet problem by graphic inte- 
gration of an integral, the integrand of which is 
a function of the Green’s function for the 
problem. Though Runge discusses the general 
theory in detail, seemingly neither he nor others 
have applied it to solution of the torsion problem. 

Though the most accurate of the graphic pro- 
cedures available for solving Dirichlet problems, 
Bairstow’s procedure is also the most elaborate 
and time consuming. Where a less accurate pro- 
cedure will suffice the entirely graphic method of 
‘‘curvilinear chequers’’ may be used to advan- 
tage. Suggested by Maxwell, it was subsequently 
developed by many writers. We mention here 
only a lengthy and exhaustive memoir by 
Richardson ;*’ the theory therein suffices for the 
solution of a great variety of two-dimensional 
problem. In a recent paper Poritsky*® illustrates 
application of this method to the torsion problem 
by effecting the stress plots of several cross 
sections. 

Perkins®® has advanced an analytic-graphic 
means of approximating the torsional moment 
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and the boundary stress. Illustratively, he ap- 
proximates M for the circle, annulus, and equi- 
lateral triangle. 


V. MISCELLANEOUS METHODS 


Of considerable practical interest is a method 
suggested by Barta® for the solution of a problem 
of equilibrium that entails determination of a 
function vanishing on the boundary of a given 
two-dimensional region. As applied to the torsion 
problem W, expressed as a Fourier-Dini series 
with unknown coefficients, is approximated by 
Y,, the sum of the first ” terms of ¥. Then suc- 
cessive substitution in ¥, of the coordinates of 
n selected points of the boundary (whereon, 
W=0); equation to zero of each of these ex- 
pressions; and solution of the resulting set of 
linear equations yield a set of preferred coef- 
ficients for ¥,. Presumably, the degree of approxi- 
imation increases with m; and proper selection 
of the points on the boundary, bilateral or quad- 
rilateral symmetry of the given domain, and 
other factors operate to afford rapid convergence 
to Vv. Thus, Barta finds that W for a square yields 
numerical values in excellent agreement with the 
exact values. 

Useful extensions of Barta’s method readily 
suggest themselves. For example, if the boundary 
of the given domain is such as to defy simple 
analytic description, it can be drafted and the 
coordinates of the selected points determined 
with an accuracy limited only by the scale of the 
drawing. Or in virtue of the shape of the given 
domain advantage may accrue to recasting the 
analysis in terms of another kind of series, say a 
Fourier-Bessel expansion. Or conformal mapping 
can be conjoined with Barta’s procedure. And 
yet other variants are feasible. 


The Fourier-Dini series has also been utilized’ 


by Csonka® in a method designed for approxi- 
mating V in a doubly connected region. Laurent’s 
theorem enables expansion of W as a sum of 
unknown harmonic functions; in turn, each of 
these functions is expressed as a Fourier-Dini 
series with unknown coefficients ; these latter are 
then determined through analytic manipulation 
of the prevailing boundary conditions; whence V 
is established. Often, as Csonka demonstrates by 
numerous examples, WV is closely approximated 
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by ¥. or V3, the sum of only the first two or three 
harmonic functions. 

A clever method (too involved to detail, here), 
invented by Schneider® for solution of the 
Dirichlet problem, conjoins use of several 
theorems of function theory. Applied to the 
torsion problem, it affords to any desired degree 
of accuracy the value of ¢ at any point of a cross 
section which can be mapped piece by piece, 
overlapping allowed, on a unit circle. Illustra- 
tively, Schneider applies his analysis to deter- 
mination of a set of values of ¢ for a cross section 
much like that treated by Runge: a cross com- 
prised of five squares, the four outermost sides 
of the cross being slightly curved. Additionally, 
he studies the ‘“‘Zweibogendreick’”’ comprised of 
the area within the outer perimeters of two 
intersecting orthogonal circles of equal radii. 

In conclusion, we mention Runge® has dis- 
proved the method advanced by Henneberg*® © 
for solution of the torsion problem. 


APPENDIX 


Though none of the following papers are concerned with 
the torsion problem, the serious student will find them of 
interest for the reasons noted: by Kryloff** and by Kan- 
torovic®? on Ritz’s method—the former contains many 
references to writings, mostly by Russian authors, not men- 
tioned in the bibliography of reference 15; by Williams,® 
Bickley,® and Frazer, Jones, and Skan7° on Galerkin’s 
method—incidentally, also on the methods of least squares 
and collocation; by Courant”! (see also part 4 of reference 
15) on the method of gradients—a variational method 
proclaimed as of considerable power, though as yet it has 
not (?) been employed for effecting actual numerical solu- 
tions; by Hencky” and by Goldsbrough"™ on aspects of a 
certain series method linked with Ritz’s method. 
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Distribution of Light Intensity within a Scattering Medium 


RoBERT R. NEWTON 
Department of Physics, University of Tennessee, Knoxville, Tennessee 
(Received June 25, 1943) 


An integral-differential equation is developed which gives the intensity within a scattering 
medium as a function of position and direction. The scattering properties of an infinitesimal 
volume of the medium are assumed known, but no special assumption is made concerning their 
nature. The equation is put in a form suitable for numerical solution. 


INTRODUCTION 


HE properties of the light scattered by a 
plane-parallel layer of a scattering medium 
have been used to indicate many of the proper- 
ties of the medium. For example, Hertel and 
Zervigon! have measured the number of fibers in 
a sample of cotton fibers which have a certain 
length by a method which essentially depends 
upon measuring the light which is transmitted by 
a layer of the sample. Sauer? has used the inten- 
sity of the light transmitted by a layer of a turbid 
medium to measure the turbidity of the layer. 
Also, many workers have used measurements of 
scattering to indicate the size of particles sus- 
pended in some other medium. *~5 
Some of the methods have depended upon 
illuminating the entire face of a layer of the 
medium, and measuring the emerging intensity, 
while, in other methods, the measurements are 
made upon the cone into which an incident 
pencil of parallel rays is spread. Only the former 
method will be discussed in this paper. 
In connection with these methods, a number 
of workers*- have developed equations which 


( 1K. L. Hertel and M. G. Zervigon, Tex. Res. 6, No. 7 
1936). 

2H. Sauer, Zeits. f. tech. Physik 12.3, 148 (1931). 

3G. I. Pokrowski, Kolloid Zeits. 47, 55 (1929). 

‘I. D. Gurewitch and G. P. Lutschinsky, Kolloid Zeits. 
60, 24 (1932). 

5G. F. A. Stutz, J. Frank. Inst. 210, 67 (1930). 

°F. Sauter, Ann. d. Physik [2] 4, 465 (1929). 

7L. Silberstein, Phil. Mag. [4] 7, 1921 (1927). 

8 A. Schuster, Astrophys. J. 21, 1 (1905). 

®Channon, Renwick, and Storr, Proc. Roy. Soc. A94, 
222 (1918). 

107. W. Ryde and B. S. Cooper, Proc. Roy. Soc. A131, 
451 (1931). 

uW. A. Fabrikant, W. L. Ginsburg, and W. L. Pulver, 
Zeits. f. Physik 81, 795 (1933). 

20. D. Chwolson, Melange Phys. et Chemie, Bull. St. 
Petersburg (1890), p. 83. 
13K. Schwarzschild, Sitz. Preuss. Akad. Wiss. (1904). 
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give the light transmitted by a plane-parallel 
layer of a medium as a function of its thickness, 
and of the optical properties of the medium. As 
a result of the complexity of the problem, sim- 
plifying assumptions have often been made. One 
assumption used is that light is scattered only in 
the direction of the incident beam, or in the 
reverse direction; another common assumption 
is that light which is once scattered is not rescat- 
tered. The first assumption will be true for 
certain types of media, e.g., a stack of glass 
plates; the second might be nearly correct for 
thin layers, but will not be met for thick layers. 

The purpose of the present paper is to obtain 
an equation which will permit calculation of the 
variation of light intensity in a layer of a scat- 
tering medium, in terms of the scattering proper- 
ties of an infinitesimal volume of the layer. In 
principle, the resulting equation can be solved, 
no matter how complicated these properties may 
be, although in practice the calculations will be 
very tedious. 

The equation obtained is similar to an equa- 
tion given in reference 11, but it is felt that in 
many cases the desired information can be 
obtained more easily from the equation given 
here. 


I, PROPERTIES OF AN INFINITESIMAL 
VOLUME OF THE MEDIUM 


Consider an infinitesimal volume of the 
medium situated at the origin of a system of 
spherical coordinates (see Fig. 1), with a thick- 
ness dz and a cross-sectional area dA. Allow light 
of intensity J to come from the direction of the 
negative z axis and strike this volume. Then the 
energy (per second) incident upon the volume 
will be IdA. Some of this energy will now suffer 
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Fic. 1. Angular relations used in defining a(@, ¢). 


‘scattering,’ using the term to indicate any 
process by which the direction of the light is 
altered. 

Assuming that the amount of energy scattered 
is proportional to the thickness of the infini- 
tesimal volume, and depends upon the direction 
which the scattered light takes, some fraction 
a(6, @)dwdz of the incident energy will be scat- 
tered into the small solid angle dw lying about 
the direction defined by (6, ¢). The energy (per 
second) thus scattered into dw will be 


IdA -a(6, )dwdz=Ia(0, ¢)dwd V, 


which depends only upon dV, and not upon the 
shape of the infinitesimal volume. The total 
energy scattered by dV will be «Jd V, where x is 
defined as: 


=f f a(@, @) sin 6déd¢. 
0 “0 


Two cases are common in practice: If the 
medium consists of small particles oriented at 
random, or spherical, the scattering will be sym- 
metrical about the direction of incidence. If the 
medium consists of cylinders with axes parallel, 
the scattered light will lie in a plane perpendicular 
to the axes of the cylinders. Considering only the 
first case, a(@, ¢) must be independent of ¢, and 
integrating once: 


c= an f a(@) sin 6d@. 
0 


It will be convenient later to use cos @ as the 
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variable, so denoting 27a(6, ¢) by o(cos @): 
c= f o(c0s 6) sin 6d8@. (1) 
0 


In addition to scattering energy, the infini- 
tesimal volume may also absorb energy of 
amount BIdV. 

In certain cases, such as for raindrops, or for 
particles whose dimensions are comparable with 
the wave-lengths of the light used, the function 
a(cos 8) could be calculated from theoretical con- 
siderations. In general, it will be easier to measure 
it by taking a thin layer of the material, illu- 
minating it by a narrow pencil of rays, and 
measuring the intensity scattered in different 
directions. That scattered directly forward or 
backward could not be measured, but could be 
obtained by interpolation. 


II. DERIVATION OF THE EQUATION 


Figure 2 represents a cross section of a plane- 
parallel layer of a scattering medium, of thick- 
ness d, and of infinite extent in the horizontal 
direction and in the direction perpendicular to 
the figure. In case the scattering particles are 
parallel cylinders, the cross section is to be taken 
perpendicular to their axes, and the extent of the 
sample along the direction of the axes is imma- 
terial. There are cases intermediate between this 
case and that in which the scattering medium 
contains randomly oriented particles, but these 
will not be considered. 

The sample is considered to be illuminated 
from above, over the entire face of the sample, 
by collimated light striking at normal incidence. 
It will be assumed that the intensity of light at 
any point within the layer is proportional to the 
incident intensity. Hence, for simplicity, the 
incident intensity is taken to be unity. 


Liititt 














Fic. 2. Cross section of a scattering layer, showing 
the coordinates used. 
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In any actual experiment, neither the sample 
nor the incident beam of light will have the 
infinite extent assumed above. A criterion for 
determining the lateral extent required will be 
given in the last section. 

The first case considered will be that in which 
the scattering particles are oriented at random. 
Referring again to Fig. 2, dA represents an 
infinitesimal area at a distance z below the 
illuminated surface. In dA lies the origin of a 
system of spherical coordinates 7, @, ¢. Further, 
dA is so oriented that its normal has the direction 
defined by the angles @ and @. 

Light which reaches a depth z consists of two 
parts: an intensity e~“+®* (per unit incident 
intensity) which has been transmitted directly, 
travelling in the forward direction, and an inten- 
sity which has reached the point in question after 
being scattered one or more times. If dV is an 
infinitesimal volume at a depth 2; below the top 
face, some of the energy reaching dA has been 
scattered from dV. 

Denote the energy of scattered light which 
strikes the area dA, and which is travelling 
within a small solid angle d2 drawn about the 
(6, ) direction, by f(z, 0,¢)dAdQ. Then f(z, 8, 
¢) denotes the energy of scattered light per unit 
area per unit solid angle which exists at a dis- 
tance z below the upper face, travelling toward 
the origin along the direction defined by (8, ¢). 

The same function at the location of dV will be 
denoted by f(2:, 41, ¢1). Then the energy which 
strikes dV with directions lying within the 
small solid angle sin @:d6:\d¢, is f(21, 01, $1) 
Xsin 6,d0,\d¢,dA,, where dA, is the effective area 
presented to light travelling in this direction. 
By the remarks of the previous section, a fraction 
(1/22)o(cos n)dwdt, of this energy will be scat- 
tered into directions such as to pass through dA. 
Here dw=dA/r’, dt; is the effective thickness of 
dV (so that dV=dAjdt,), and 7 is the angle 
between the directions defined by (6, ¢) and by 
(91, @1), so that 


cos n=cos 6 cos 6;+sin @ sin 6; cos (6—¢1). 


The energy scattered from the small solid angle 
sin 6,d0\d¢; into dw is then (1/22)e(cos n) 
X f(21, 01, 61) sin 0:d0,d¢id Vdw. Upon integrating 
this over 6; and ¢;, there is obtained the total 
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energy which is rescattered by d V into the proper 
direction to pass through dA. 

In addition, the intensity exp [—(«+ 8)2: ] 
reaches dV without being previously scattered ; 
of this, an amount 


exp [ —(xk+8)21 ][o(cos 6)/2x |d Vdw 


is scattered so as to pass through dA. Then the 
total energy scattered by dV into the solid angle 
dw is: 


(1/ 2) f f ofcns @ifite: Cx. ds) ini Olid 
0 0 
+exp (—-yz1)0(cos 02 Vdw, 


writing y for (x+ 8). Of this energy, a fraction 
e-7" will arrive at dA; the rest will be absorbed 
or scattered out of the path. That which reaches 
dA will be travelling in directions lying within 
dQ, the solid angle subtended by dV at dA. This 
energy gives a contribution to f(z, 6, ¢)dAdQ, 
which we call df(z, 0, ¢)dAdQ. Taking dV to be 
a pillbox with axis parallel to 7, of thickness dr, 
then dV =r°dQdr, while dA =r°dw. Making all of 
the indicated substitutions, we obtain the 
equation : 


df(z, 0, o)r*dQdw 


= a/2=)| f f a(cos ) f(21, 81, d1) sin 0:d0,dq, 
0 0 
+exp (—y21)¢(cos 0) Jerrad (2) 


This can be divided through by r°dQdw, and 
integrated over r from zero to the boundary of 
the layer. Further, the function f(z, 0, ¢) must 
be independent of ¢. This gives the integral 
equation for the function f(z, 6): 


fle, 0) = (1/2) f f F f “e(cos 9) flu 01) 


Xsin 0:d0,d¢,+exp ( — y21)¢(cos 0) fewer (3) 


The various ranges for r are: for 0£< @<7/2, 
0<rXz/cos 6; for 0=2/2, OS r& «; for r/2<@ 
<r, OX r K (d—2z)/cos 0. 
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Changing the variable to z; by means of the 
relation r= (z—2,)/cos @: 


f(z, 0) = (exp [ —yz/cos 6] /2x cos 8) 


«ff f a(cos 7) f(z1, 41) 


Xexp [(yz:/cos 6) ] sin 0:d0:d¢, 


+ (cos 6) exp [y2:(1—cos @)/cos 6] as (4) 
for 0642/2, while for 6= 2/2: 


f(z, 2)=(1/2e0)| f f a(cos ») f(z, 81) 
sin Odbudon+o(0)e-™| (5) 


The lower limit of integration in the first equa- 
tion is zero, if 0<2/2, or d, if @>-7/2. 

By differentiating Eq. (4) with respect to z, the 
result can be put in a simpler form, and one 
which also holds for 6= 2/2, so that the trouble- 
some double form is eliminated. Performing the 
differentiation and solving for f(z, 6) : 


f(z, 0) = —(cos 0 ‘y) [Of (sz, 0)/dz] 
+[a(cos 0)/2ry je~” 


+(1, 2m) [ f a(cos n) f(z, 41) sin 0,d0\dq,. (6) 
0 0 


The problem for the case in which the scatter- 
ing medium consists of small cylinders with axes 
parallel is handled in much the same manner, 
and only the result will be given. The problem 
is only two-dimensional, so that o(cos @)d@dz is 
the fraction of the energy striking it which a 
small thickness dz of the medium scatters into 
a wedge with dihedral angle dé. The figure 
corresponding to Fig. 2 is a cross section of the 
medium taken perpendicular to the axes of the 
cylinders, so that all directions in which energy 
is travelling lie in this plane. If 


=f a(cos 0)d@ 


denotes the scattering coefficient, and 8 the 
absorption coefficient, the equations correspond- 
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ing to Eqs. (3) and (6), respectively, are: 


f(z, 6) -flf o(cos n) f (21, 6:)d0; 
+exp (—y21)¢(cos 0) evar, (7) 


f(z, 0) = — (cos 6/y) (f(z, 6) /dz] 


+[o(cos 6) /y je-” 
+(1/y) f o(cos n) f(z, 0:)d0;, (8) 


where cos n=cos (#@—6,); the definitions of 
other quantities, and the ranges of integration 
for r, are as before. 


Ill. USE OF “FORWARD” AND 
“BACKWARD” FLUX 


Several authors’~'® have treated the problem 
by assuming that all of the scattered energy 
travelling in directions for which @ lies between 
0 and 2/2 is equivalent to a forward flux F,, 
travelling in the direction of the incident light, 
and that the energy with directions for which 6 
is greater than 2/2 is equivalent to a backward 
flux F2. Then light is assumed to be scattered 
either forward with a coefficient ¢x, or backward 
with a coefficient (1—¢)«. Ryde and Cooper'® 
have elaborated this treatment further by using 
different coefficients ¢ and ¢’ for the parallel, 
directly transmitted, light, and for the scattered, 
diffused, light. In deriving the coefficient for the 
diffused light, however, they assumed the diffused 
light to be distributed equally in direction, 
whereas the coefficient would actually have to 
depend upon the angular distribution of the 
light, and hence upon position within the scat- 
tering layer. 

Equation (8) can be reduced to the form given 
by Silberstein,’ if o(cos @) and f(z, @) have the 
following forms: 


For —7/2<06< 2/2, a(cos 0) =£xd(8), 
f(z, 0) = F,6(6). 
For 2/2<0@<3nr/2, o(cos 6)=(1—f)xi(@—7Z), 


f(z, 0) = F.5(0—7). 
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For convenience in making these definitions, the 
range of @ has been changed from (—z7, 7) to 
(—2/2, 3/2). 5(0@) is the delta function,'4 which 
has the properties: 


i(x—a)=0, x¥a, 
and 


f see —a)dx = f(a), 


provided that the range of integration includes 
the point a. Substituting the values of o and 
f(z, 0), we obtain, for || <r/2: 


+ F\8(8) = —cos 0Fy'5(0) + ¢xe-*75(8) 


+ f ¢x8(0—61) F,8(0:)d0, 


3n/2 
+f (1 —£)x6(0— 0, — 2) F25(0; — 3) dy. 


Evaluating the integrals: 
+ F\6(8) = —cos 0Fy'5(6) + ¢xe-7"6(8) 
+ £«F16(0) + (1—¢) x F26(6). 
Integrating over 6 from —2/2 to 7/2, 
yFy= — Fi +fcne-V + o«KFit+(1—f)K Fe. (9a) 
Proceeding similarly for the case 7/2 £6< 3x/2, 
yP.= Fo +(1—f$)xe-% +(1—$) x FitgnFe. (9b) 


If Eqs. (9a) and (9b) are solved for the deriva- 
tives F;’ and F,’, they are seen to be identical 
with the differential equations given by Silber- 
stein, of which the integration is given in 
reference 7. | 


IV. DISCUSSION 


The appearance of the term (¢(cos 6) /2mry)e~”* 
in Eq. (6) is interesting, for if ¢ were a constant 
equal to «/2, this term becomes (1/47)e~™, 
neglecting the absorption coefficient compared to 
the scattering coefficient. This represents a total 
intensity at any point of e~“, distributed evenly 
in direction ; and represents a first approximation 
to the solution of the problem. Similar remarks 
apply to Eq. (8). 





“For a discussion of this function see P. A. M. Dirac, 
Quantum Mechanics (University Press, Oxford, 1935), sec- 
ond edition, Chapter IV. 
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The case for which o equals «/2, and there is 
no absorption, so that y=x, is interesting from 
another consideration. Upon making these sub- 
stitutions in Eq. (6), it can readily be seen that 
the equation contains only the combination xz; 
that is, xz is a dimensionless variable, and xd is 
a dimensionless parameter of the system which 
determines how much light is transmitted by the 
medium without suffering scattering. Under 
other conditions, y always appears as yz, or as 
(o/y). If the total number of particles in a scat- 
tering layer remains unchanged, these quantities 
remain unchanged if the layer is compressed or 
expanded. In this sense, yd is a parameter of the 
system, with the same meaning as xd has above. 

The form of Eqs. (6) and (8) is such that they 
will probably defy analytic solution, even for 
very simple forms of the function o(cos @). It will 
then be necessary to obtain the function f(z, @) 
by methods of successive numerical approxima- 
tion. If absorption can be neglected, the equa- 
tions are somewhat simplified ; and if, in addition, 
the scattering can be assumed spherically sym- 
metrical, the resulting equations are consider- 
ably simpler. The integral 


f a(cos n)d¢ 
0 


can be evaluated as a function of @ and 6; as soon 
as the form of o is known, so that only one 
integration need be performed for each step in 
the approximation process. In cases where it 
is felt that the derivative (6f/5z) cannot be ob- 
tained numerically with sufficient accuracy, the 
double integral form of Eq. (4) can be used. 

From Eq. (3), an approximate criterion for the 
lateral extent required of the sample can be 
developed. Referring to Eq. (3), the quantity in 
square brackets may be considered to behave 
roughly like exp( —yz:). Then, if a designates the 
horizontal distance from the origin of r to the 
edge of the layer, the error in calculating f(z, @) 
will be less than one percent if: 


z/cos 6 
a | e777 (l—cos dr 
0 


z/cos 6 
2 100-e-% f grits ey, 
a 


/sin 6 
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since 23+r=z2+r(1—cos 6). That is, the con- 
tribution to the integral from points beyond the 
actual edge of the layer is less than one percent 
of the entire integral. Evaluating the integrals, 
it may be seen that the largest value of a required 
is for 6=90°. Substituting this value of @, the 
above condition becomes: 


e~7* £0.01. 


The necessary lateral extent is then independent 
of the thickness of the sample. Of course, actually 
cutting away the portion of the sample lying at 
a greater distance than a from the point where 
measurements are being made does not produce 
the same effect as ignoring that part in the cal- 
culations. In spite of this, the criterion given 
should prove useful. 


There are several uses for the function f(z, 6) 
which suggest themselves. If one is interested in 
the directed intensity at any point, that is, of 
course, just f(z, 6) itself. The emerging intensity 
in any direction is the value of f(z, 0) at z=d. 
To obtain the total intensity received upon a 
flat surface, it will be necessary to multiply 
f(z, 9) by the cosine of the appropriate angle, and 
to integrate over all directions from which light 
is received. This process, if carried out for a 
surface coincident with the lower face of the 
sample, and e~*4 added to the result, will give the 
total transparency of the layer. 

The author wishes to express his appreciation 
to Dr. K. L. Hertel, of this department, who sug- 
gested the problem, and to Dr. W. G. Pollard, 
who has contributed many valuable suggestions. 
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The differential equations governing the deflection of structural elements such as loaded 
cords and beams freely suspended or resting on an elastic foundation are solved operationally. 
The operational method simplifies the analysis because it reduces the number of arbitrary 
constants to be determined and enables the effect of concentrated loads or applied couples to be 


determined directly. 


INTRODUCTION 


HE determination of the equilibrium con- 

figurations of loaded cords and beams leads 
to the solution of ordinary differential equations 
subject to multi-point boundary conditions. In 
the usual applications, however, the boundary 
conditions are specified at two points. Many 
analytical and graphical methods! have been 
devised to solve problems of this nature. In 
some instances, purely analytical procedures are 
very cumbersome. It is the purpose of this dis- 
cussion to illustrate how problems of this type 
may be solved simply and directly by the use 


1S. Timoshenko, Strength of Materials (D. Van Nostrand 
Company, New York, 1932). 
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of the Laplace transform or operational method. 
It may be mentioned that in the usual applica- 
tions of this method, the differential equations 
under consideration are solved subject to one- 
point boundary conditions; however, the utility 
of this method in certain multi-point boundary 
value problems will be apparent from this dis- 
cussion. 


I. THE LAPLACIAN TRANSFORMATION 
a. General Considerations 


Before proceeding, it may be well to review 
briefly the theory on which the modern opera- 
tional calculus is based.2 The modern theory of 


2L. A. Pipes, J. App. Phys. 10 (1939). 
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the operational calculus, a mathematical tech- 
nique which has proved to be such a powerful 
mathematical tool in the solution of the linear 
differential equations of applied mathematics is 
based on the Laplacian transformation. This 
transformation is usually written in the following 
form: 


g(p) =p e ’*h(x)dx Re p>0o. (1) 


In the usual applications of the theory, A(x) is 
a function for which the definite integral (1) 
exists and satisfies the condition, 


h(x)=0 for x<0. (2) 


A rigorous discussion of the restrictions that 
must be imposed on a function h(x) in order that 
it may have a Laplacian transformation is given 
elsewhere. It may be said, however, that the 
functions encountered in the various physical 
problems where the transform method is applica- 
ble satisfy the required conditions. In most of the 
modern literature on operational or Laplacian 
transform methods, the functional relation ex- 
pressed between g(p) and h(x) is written in the 
following form: 


g(p) =Lh(x). (3) 


The Z denotes the ‘Laplacian transform of”’ 
and greatly shortens the writing. The relation 
between h(x) and g(p) is also written in the form 


h(x) =L~'g(p). (4) 


In such a case we say that h(x) is the inverse 
Laplacian transform of h(x). 

A table of definite integrals of the form (1) 
may be used to obtain many transforms. A col- 
lection of these transforms useful in applications 
to the theory of structures will be found in the 
appendix. 


b. The Laplacian Transforms of Derivatives 


In the solution of differential equations by the 
Laplacian transformation, it is necessary to 
transform the derivatives of a function. It is 


3N.W. McLachlan, Complex Variable and Operational 
Calculus (Cambridge, 1939). 

*M. F. Gardner and J. L. Barnes, Transients in Linear 
Systems (Wiley, 1942). 
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easy to show by integrating by parts? that the 
following relations exist: 


If 


Ly(x)= Y(p), (5) 
then 


L(dy/dx)=pY—pyo, (6) 
L(@y/dx*) =p? Y — p*yo— py, (7) 
L(d*y/dx*) = p* Y — p*yo— p*y1— pyr, (8) 
L(d*y/dx*) = p'Y — p'yo— p*yi— P*y2— pys, (9) 


where 


y,=(d'y/dx’) evaluated at x=0. (10) 


These relations are of great importance in the 
solution of differential equations by the use of 
the Laplacian transformation. 


II. DEFLECTION: OF LOADED FLEXIBLE CORDS 


Let us consider the deflection of the tightly 
stretched cord of Fig. 1. 


waar (2) 
LJ A444 141147 4 » 








Fic. 1. 


It may be shown? that the equation governing 
the small deflection y of a cord tightly stretched 
by a tension T and loaded by a load of w(x) per 
unit length is 


d*y/dx*? = —w(x)/T, (11) 


where y is measured downward and w(x) is 
positive when acting in a downward direction. 
To solve this equation, let us introduce 


Ly(x)=Y(p), (12) 
Lw(x) =g(p). (13) 
Then by using (7), Eq. (11) is transformed to 


Pp? Y — p*yo— py = —g(p)/T. (14) 


1. Uniform Load 


To illustrate the method, let us first consider 
the case in which 


w(x) =Wo. (15) 


5 T. Von Karman and M. A. Biot, Mathematical Methods 
in Engineering (McGraw-Hill, 1940), Chapter 7. 
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In this case (see Appendix), we have 
Lwy = Wo. ( 16) 


If the two ends of the span are at equal heights, 
then 

y=Oatx=0, y=Oatx=s. (17) 
Hence 


yo=0. (18) 
Equation (14) then becomes 
Y=(y1/P) —(wo/T)(1/p*). (19) 


By the table of transforms in the appendix, the 
inverse transform of (15) is 


Y= ix — (Wo/T)(x*/2). (20) 


To determine y; we make use of the boundary 
Eq. (17) and obtain 


V= (Wo ‘T) (s, 2). (21) 
Substituting this into (16) we obtain 
y = (wo /2T) (sx —x?) (22) 


for the deflection of the cord. 


2. Uniform Load Extending over Part of Span 


Let us consider the case of loading shown in 
Fig. 2. In this case, the transform of the load is 





K=O ao) 
y y 
%, a, V 
bs 
Fic. 2. 
given by (see Appendix) 
g(p) = wo(e~ 7” — €72”), (23) 
Equation (14) becomes 
>, wir <-¢Cc 
9 ee ee -), (24) 
p T p* 


The inverse transform of this is (see Appendix) 


y=yix O<x<Kx1, 
y=Vix—(wWo/2T)(x—x1)?  x1<K<x2, (25) 
Y= ix — (Wo/2T)(x—-x1)? 

+(wo/2T)(x—x2)? xe<x<s. 
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To determine y; we make use of the boundary 
condition 


y=0 at x=s. (26) 


Inserting this in the third equation of (25) we 
obtain 


y= (wo 2Ts) [(s —x;)? oie (s —X2)* ]. (27) 
The deflection curve is given by (25) with this 
value of yi. 
3. Concentrated Load 


Let us compute the deflection of a flexible cord 
fixed at both ends and supporting a concen- 
trafed load W at a distance of x from the left 
end as shown in Fig. 3. 


hg 
Z=0 Aa=S 
y | y 

A ~, V 


Fic. 3. 





To determine this deflection, we may use Eq. 
(14). It is shown in the appendix that the trans- 
form of a concentrated load W situated at x=x, 
is given by 


g(p) = We-*'”p. (28) 
Substituting this into (14), we obtain 


1 W 1 
¥ =-y,;——e ™?-, (29) 
p 1 p 
Computing the inverse transform by using the 
table in the appendix, we obtain, 
y=yix O0<x<x, (30) 
y=yix—(W/T)(x—x1) x1<x<s. 


The constant y; may be determined from (30) by 
realizing that 


y=0 at x=s. (31) 
Hence we obtain 
yi=(W/Ts)(s—%x)). (32) 
Substituting this value of y; into (30) we ob- 
tain the following function for the deflection. 


y=(W/T)x[1—(x1/s)] O<x<x, (33) 
y=(W/T)xif1—(«/s)] x1<x<s. 
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4. Arbitrary Load 


To determine the deflection produced by an 
arbitrary load w(x) per unit length, w(x) we 
have from Eq. (14) 


1 _1 gd) 


=—Y) " 


p T Pp 





(34) 


To compute the inverse transform of g(p)/p* 
we use the Faltung theorem (see Appendix) ; we 
thus obtain 


g(p) 7 
n= f w(u)(x—u)du. (35) 
p° 0 
The inverse transform of (34) is 
1 z 
iat a f w(u)(x—u)du. (36) 
0 


The condition that y must vanish at x=s 
enables y; to be determined. Its value is 


1 8 
y=— f w(u)(s—u)du. (37) 
Ts 0 


Substituting this value of y; into (36) we obtain 


8 


x 
y= J w(s—wdu 
Ts 0 
1 z 
<= f w(u)(e— udu (38) 


0 


for the deflection of the cord under the influence 
of a general load w(x). 


III. STRETCHED CORD WITH ELASTIC SUPPORT 


Let us consider the deflection of a cord stretched 
by means of a tension JT resting on an elastic 
foundation and under the influence of a certain 
load w(x) per unit length. 

The differential equation governing the deflec- 
tion of such a string is given by® 


€ 


T——ky=—w(x). 


(39) 
dx? 


In this equation, T is the tension of the cord, 
w(x) is the load per unit length, and k a constant 
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of proportionality which depends on the stiffness 
of the assumed elastic foundation. 

That is, the vertical deflection of the cord is 
restrained by a large number of springs such that 
their action can be replaced by a distributed 
restoring force per unit length equal to ky. 

Equation (39) may be written in the form 








d*y — w(x) 
pe —Ccy= , (40) 
where 
e=k/T. (41) 
To solve this equation operationally, let 
Ly(x)=Y(p), (42) 
Lw(x)=g(p). (43) 
Equation (40) is thus transformed to 
(p?—C)V=pyotpyi—Lg(b)/T]. (44) 


1. Cord Fixed at Its End, Concentrated Load 


Let us consider the case shown in Fig. 4. 


‘X, 


—— 














, i hk 


1=0 A=s 


Fic. 4. 


In this case, we wish to determine the vertical 
deflection y produced by the application of a 
concentrated load W at a distance x; from the 
left end support. The cord is under a tension T 
and rests on an elastic foundation of elastic 
constant k. 

We have the following boundary conditions: 








yo20, y=0 at x=s. (45) 
The transform of the concentrated load at 
x=, is 
g(p) =e *'"pW. (46) 
Hence Eq. (44) becomes 
py W 
¥o—_—_-_¢~——, ij 
(p?—c*) 7 (p? —c*) 
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Consulting the table of transforms in the 
appendix, we find the inverse transform to be 


sinh (cx) 


y=" - O<x<x, (48) 
c 
sinh (cx) W 
y=" —— sinh [c(x—~x,) ] 
c Tc 


x,<x<s. (49) 


If we use the boundary condition, that the 
deflection must vanish at x=s, we obtain 
W sinh [c(s—x;) ] 
n° ——_————_——. (50) 
7 sinh (cs) 
Substituting this value of y; in Eqs. (48) and (49), 
we obtain the following expression for the deflec- 
tion of the cord. 


W sinh [¢(s—x,) ] sinh (cx) 


y= : ". O<x <x, 
Tc sinh (cs) 





(51) 
W sinh [e(s—-x1) ] 


— sinh (cx) 








y = 
Tc sinh (cs) 


——-sinh [c(x—x1)] x1<x<s. 
Tc 


2. Infinite Cord Elastically Supported, 
Concentrated Load 


The deflection of an infinitely long elastically 
supported cord under the influence of a concen- 
trated load W may be obtained as a special case 
of (41). 

In order to obtain the deflection in this case, 
let 

z=(x—%), (52) 


x,=s5/2. (53) 


Making these substitutions in (51), we obtain 


W cs cs W 
y =— sinh ( ) sinh ( +c) sinh (cz), 
Tc 2 2 Tc 


z>0. (54) 
Then 
— W W 
Limit y=——e-“ =e" _ z>0. (55) 
— 2Tc 2(kt)! 
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This is the required deflection and z is meas- 
ured from the point of application of the load. 


3. Elastically Supported Cord with Uniform Load 


Let us consider the case shown in Fig. 5. 


oe a 


/ 
V7RKTTTi777TI7V 
=o ' ' ass 
<4 a 


‘ a 





=z 
Fic. 5. 


In this case, the cord is loaded with a uniform 
load wo per unit length. The load extends from 
X1 to Xo. 

In this case, we have 


Lw(x) = Wo(e-71? —€ 72?) , (56) 
Equation (44) becomes 


p Wo f € 7? — e- 72P 
Y=yi— —-—( ewes | ). (57) 
(pp—c?) 7 





Consulting the table of transforms in the ap- 
pendix, we obtain 


sinh (cx) 
y=Vir— —— O<x<x1, 
c 


sinh (cx) Wo 
y= -———f{cosh c(x—x1)—1} 
c Tc? 
X1<X<Ne, (58) 


sinh (cx) Wo 
y a y1 "ae 7 9 
c ic 


cosh ¢c(x— 1) 


Wo 
+—coshc(x—x2) x2<x<s. 
Tc? 


Using the boundary condition that the deflec- 
tion must vanish at x =s, we obtain the following 
value for the constant y;: 


Wo 
n= {cosh c(s — x1) 


Tc sinh (cs) 





—cosh c(s—x2)}. (59) 


Substituting this value of y; into (58) gives the 
required deflection. 
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4. Elastically Supported Cord, General Loading 


The transform of the deflection of an elastically 
supported cord under the influence of a general 
load w(x), is by Eq. (44) 








1 
y= Py 1 a) 7 (60) 
(p?—cp) T (p?—c*) 
where 
Lw(x) =g(p). (61) 


By the Faltung theorem in the appendix we have 


1 z 
+= =f w(u) sinh [c(x—u) du. (62) 





4 


The inverse transform of Y is thus given by 


sinh (cx) 1 # 


y=yy-———_ = f w(u) sinh [c(x—) du. 


‘ 


Cc C#Ho (63) 


The fact that the deflection vanishes at x=s 
leads to the following value of -y;: 





Vi= 


1 8 
~ f ww) sinh [c(s—u) ]du. (64) 


T sinh (cs) 
Substituting this value of y; into (63) gives the 


required deflection. 


IV. THE DEFLECTION OF BEAMS BY 
TRANSVERSE FORCES 


The deflection y of a uniform beam of flexural 
rigidity EJ satisfies the following differential 
equations :! 5 


EI (d*y/dx*) = w(x) =load per unit length. (65) 
F= — EI (d*y/dx*) =the shear force. (66) 
M=-—EI(d*y/dx*) =the bending moment. (67) 


If the load per unit length is known, then the 
deflection y may be determined by solving Eq. 
(65) subject to the boundary conditions of the 
particular problem under consideration. In order 
to solve (65) operationally, let 


Ly(x) = Y(p), (68) 
Lw(x) = g(p). (69) 


Using Eq. (9) for the transform of the fourth 
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derivative, Eq. (65) may be transformed into 


- g(p) 
EI p* 





Yi yo Y3 
Y a ca awe (70) 
db Pp p® 
This equation is basic in the calculation of the 
deflection of beams. 


1. Uniform Beam Clamped Horizontally at Both 
Ends Under the Influence of a Uniformly 
Distributed Load 


Consider the deflection of the beam shown in 
Fig. 6. 














q 
( 
X=S 


Fic. 6. 


In this case, the deflection and slope at x=0 
are both equal to zero, hence 


yor y=. (71) 
The transform of the uniform load wo is given by 
g(p) =Wo. (72) 
Hence Eq. (70) becomes 
Wo 1 yo Ys 
Y=— —+—+— (73) 
EI p* Pp? p® 
and 
Wo x* x ro 
ee —+¥e- +%s—. (74) 
EI 4! 2! 3! 


The constants y2 and ys are determined by the 
condition 


y=dy/dx=0 atx=s. (75) 
These conditions yield 

yi= —(1/12)(wos/ED). (76) 

yo= (1/12) (wos?/ ED). (77) 


Substituting these values of y; and ye into (74) 
we obtain the following equation for the de- 
flection : 

y = (1/24) (wo/ EI) x?(s —x)?. (78) 


2. Uniform Beam Clamped Horizontally at Both 
Ends and Carrying a Concentrated Load 


Consider the problem of determining the de- 
flection of the beam shown in Fig. 7. 


491 




















In this case, the transform of the concentrated 
load is 


g(p) = Wipe”. (79) 


Since the deflection and slope are zero at 
x=s, Eq. (70) for this case becomes 


W €~ 71P v2 ¥3 


¥=——4—4— (80) 
EI p* Pp p' 
The inverse transform of Y is 
=* x* 
V=Vo—+Y¥3s— O0<x <2, 
2! 3! 
(81) 
W x* x 
y=—(x—x1)? + y2—+y3s— x11<x<s. 
6EI 2 6 


The conditions that the deflection and slope 
must vanish at x=s enables y, and y2 to be de- 
termined. Inserting these values of ys and y; into 
(81) we obtain the following equations for the 


deflection : 
1 W x? 
y =—- — —(s—x1)*{3xys—(s+2x1)x} O<x<m, 
6 EI s' 
(82) 
1 W x" 
y =—- — —(s—x)*{(3s—2x,)x—xis}  xi<x<s. 
6 EI s* 


3. Uniform Beam Clamped Horizontally at One 
End and Free at the Other Carrying a 
Concentrated Load 


- As another example of the general method, 
consider the problem of determining the deflec- 
tion of the beam shown in Fig. 8. 




















a 
W 
7 
cee > 
/ 
Fic. 8 
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In this case, the basic Eq. (70) becomes 


WwW €7 oP V2 ¥3 


Y=— ——_+—+~—. (83) 
EI p® pP p® 
The inverse transform of Y is 
X2 x? 
Y=Vo—+¥3— O<x<a, 
2 6 
(84) 
W x? x* 
y=———-(x—a)? + yo—+33--  a<x<s. 
6EI 2 6 


In this case, since the end x=s is free, it 
follows that the bending moment and shear 
force at that point must vanish. By (66) and 
(67), we therefore have 


d*y/dx?=d*y/dx*=0 atx=s. (85) 
Inserting these conditions in (84), we obtain 

y2=aW/EI!, (86) 

ys= —W/EI. (87) 


Substituting these values into (84) we obtain 
the following equation for the deflection : 


Wx*/a x 
y-——(“-") O0<x<a, 


EI \2 6 
(88) 
Wa*/x a 
y-—(=-*) a<cx<s. 
EI \2 6 


4. Uniform Beam Clamped Horizontally at Both 
Ends Under the Influence of an Applied Couple 


Consider the case shown in Fig. 9. 














S 


XK=S 


Fic. 9. 


It is required to determine the deflection of 
the beam produced by the application of a 
couple My at x=a. The application of the couple 
may be regarded as the limit of the loading 
shown in Fig. 10. 
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a mi 


I, 


Fic. 10. 














i. 


That is, the couple may be regarded as the 
limiting form of loading as shown in Fig. 10 
provided that 

Limit (F6) = Mo. 


Foo 
6-0 


(89) 


To calculate the deflection in this case, we 
write Eq. (70) in the form: 


Y=— —_+—+ (90) 


where g(p) is the transform of the loading shown 
in Fig. 10, that is 


g(p) = — pFe-*?+ pFe““t)P (91) 
= pFe—*?(e—§?—1). 
Now 
Limit g(p) = Limit — (F8)p*e-*”, (92) 
6—0 5-0 
F-2 Fox 
= — Mop*e°?. 


This is the transform of a couple of strength Mo 
impressed at x=a. Substituting this value of 
g(p) in Eq. (90), we obtain 


Yo ———+4—+—., (93) 
EI p? p* p 
The inverse transform of y; is 
x? -y3x3 
vY=Vo—-+ O0<x<a, 
z.6hCUG 
(94) 
— My (x—a)? x* x* 
y = — —FH+ 2-H 93— a<xc<s. 
EI 2 2 6 


The condition that the slope and the deflection 
must vanish at x=s leads to the following values 
of the constants y2 and ys: 











Mo 
yo =——{s?—4as+3a’}, (95) 
SEI 
6M,a(s —a) 
y= cee (96) 
SEI 
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Substituting these values of yz and y3 in Eq. 
(94) gives the required deflection. 


V. DEFLECTION OF BEAMS ON 
ELASTIC FOUNDATION 


Let us consider the deflection of a uniform 
beam that is attached to a fixed base by means of 
an elastic medium. The deflection y of a beam of 
flexural rigidity EJ resting on such a foundation 
satisfies the differential equation! 5 


EI (d*y/dx*) +ky=w(x). (97) 


The factor k is called the modulus of the 
foundation and depends on the elastic property 
of the foundation. w(x) is the load per unit 
length applied to the beam. For simplicity, let 
us write 


a=(k/4ED)}. (98) 


Equation (97) may then be written in the form 
(d*y/dx*) +4aty =w(x)/EI. (99) 


To solve this equation operationally, we write 


Y=Ly(x), (100) 
g(p) = Lw(x). (101) 
Equation (99) is then transformed to 
(p'+4a*) Y=[¢(p)/EI]+p*yo 
+P tp y2t+pys. (102) 


This equation is basic in the computation of 
the deflection of loaded beams on elastic foun- 
dations. 


1. Beam on Elastic Foundation Clamped Hori- 
zontally at Both Ends Under the Influence 
of a Concentrated Load 


Let us consider the deflection of the beam 
shown in Fig. 11. 


A |W 


ip 








~~ 








,TTTTTTTTTTTTTT7 








/ 
azo a=s 
Fic. 11. 
In this case, the conditions that the slope and 
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deflection at x=0 must vanish lead to 
yo= yi=0. (103) 


The transform of the concentrated load W, 
acting at x=), is 


g(p) = Wpe-”®. (104) 
Hence in this case, Eq. (102) becomes 
W pe p? p 
TE spe. oF . pene 
EI (p'+4a*) (p! +4a! ) (p*+4a‘) 
(105) 
= Ly. 


The transform of dy/dx is given by 
dy Wp” p* 

* 

dx/ EI(p*+4a*) (p*+4a*) 











— . (106) 
* (pt4-4a") 
Now let 
p 
L- ——=9¢1(x), (107) 
(p*+4a*) 
pet, (108) 
(p'+4a‘) 
3 
| a — = 3(x), (109) 
(p*+4a‘) 
4 
L“ : —= 4(x). (110) 





(p*+4a‘) 


The values of ¢,(x) are found in the Appendix. 
In terms of these functions, we have the inverse 
of Y given by (105) in the form 


Y=YVoh2(x)+ysb1(x) O<x<b, 
y= (W/EI)$1(x — b) + y2g2(x) +y361(x) 
b<x<s. 


(111) 


The constants y, and ys may be found from 
the condition that the deflection and the slope 
must vanish at x=s. From (111) we thus obtain 


= (W/EI)¢1(s —b) +y2o2(s) +yadi(s). (112) 
And from the inverse of (106) we have 
0= (W/EI)$2(s—b) +y2os(s) +yab2(s). (113) 
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These two equations may be solved for the 
constants ye and y3. This gives 


W [o1(s)o2 (s—b) —¢: $2(s)1(s — 5) ] 


y2= — (114) 
EI [-$2°(s) — $a(s)du(s) ] 
W [$2(s)$2(s —b) — o3(s)oi(s — 6) ] 

V3 =  ——_..._ (11) 
Er [o1(s)bs(s) — oe 2(s)] 

The deflection is obtained by substituting 


these values of yo and ys; into (111). 


2. Beam on Elastic Foundation, Clamped Hori- 
zontally at Both Ends Under the Influence 
of a General Load 


In this case, since the deflection and slope 
vanish at x=0, Eq. (102) becomes 

















1 g(p) p* ye rae PYs 
=—_ = Vv. 
EI (p'+-4a") (p'+4a*) (p'+4a*) (116) 
The inverse transform of the slope is given by 
1(2)- 1 pg() 
dx EI (p'+4a') 
3Vo =; 
5 
(p'+4a*) (p*+4a*) 
By the Faltung theorem, we have 
g(p) ¥ 
yt = [ wiudor(e—wu, (118) 
(p*+4a*) 0 
p z 
L-\_—__= f w(w)oa(x—u)du. (119) 
(p'+4a*) 0 


From (116) and (117) we thus obtain the de- 
flection and slope to be 


1 z 
y=— f w(u)doi(x—u)du 
EI Jo 


+yob2(x) +yabi(x), (120) 
d 1 z 
ep J Ue (x—u)du 

+ yohs(x) +¥32(x). (121) 
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It we let 








1 8 

A= J wou(s—wdu, (122) 
EI Jo 
1 8 

B= J w(woa(s—u)du, (123) 
EI 46 


and make use of the fact that the slope and de- 
flection both vanish at x=s, we obtain the 
following values for the constants y2 and ys: 


Boi(s) —A¢2(s) 
$2°(s) —o1(s)$a(s) 
7 Bo2(s) —A¢a(s) 
+e G1(s)oa(8) — 62°%(5) 





y= 


(124) 





(125) 


If these values of ye and y3 are substituted 
into (120), the value of the deflection is obtained. 


VI. GENERAL CONSIDERATIONS 


The above examples show the utility of the 
operational or Laplacian transform method of 
solving the differential equations which occur in 
the theory of structures. 

The examples given are typical and other 
types of boundary conditions may be solved in 
the same manner. The advantages of the opera- 
tional method are that one may write one equa- 
tion for the whole span and. also the work of 
determining the arbitrary constants of the solu- 
tion is considerably reduced. 


APPENDIX A. BASIC THEOREMS OF THE 
LAPLACIAN TRANSFORMATION 


The utility of the Laplacian transformation in the solu- 
tion of differential equations is based on some important 
relations which follow as a consequence of the fundamental 
Eq. (1). The most important of these theorems will be listed 
here for reference. They are established in references 2-4. 


If Lh(x) =g(p) then 


Lh(sx) =g(p/s) I. 
where s is a constant s>0. 
L(dh/dx) = pg(p) — ph(0), Il. 
where h(0) is the value of h(t) at t=0. 
d"h k=n-l 
(~~) =pre) — > AhOpr, IIL 
dx” k=0 
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where h!*)(0) denotes the value of the kth derivative of 
h(x) evaluated at x=0. 


(0 x<a 
L~'e~4?g(p) "7 . IV. 
h(x—a) x>0 
where a is a positive constant. 
If L-g:(p) =hi(x) and Lo'g2(p) =h2(x), then 
£:(p)g2(P) edie, 
L-1\—__—__ = V. 


p | fozhe(u)hi(x—u)du 
This theorem is known as the Faltung theorem. 
APPENDIX B. TABLE OF LAPLACE TRANSFORMS 


g(p) =p [<P x) dx 

















No. g(p) h(x) 
1 p 5(x), the Dirac impulse function 
5(x) = © for x=0, 5(x) =0, x ~0 
[sede =| 
1 x” gt le 
2. p ni, 2,3, --> “ (n a positive integer) 
3 on cos (ax) 
4, Ces sin (ax) 
8 aon cosh (ax) 
6. gai sinh (ax) 
1 1 
a (pP—a?) a2 cosh (ax) —1 } 
a 
8. ortas qqacsin ax cosh ax 
—cos ax sinh ax) =¢;(x) 
9. ort pa(sin ax sinh ax) =¢2(x) 
= . dig 4 
10. (p'+-4a") 34 sin ax cosh ax 
+cos ax sinh ax) = ¢3(x) 
4 
e2. = 2s cos ax cosh ax = d(x) 
(p'+4a‘) 
aie Tie apie ake 
12. (p—a 2q3¢sinh ax—sin ax) 
- a — 
13. Ga) rs osh ax—cos ax) 
14. eon | tials ax+sin ax) 
(p'—a’*) 2a 
15. - - *(cosh ax+cos ax) 


(p'—a*) 2 
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Radioactive mercury, Hg'*?, was prepared by deuteron bombardment of gold. This isotope 
was used as a tracer to establish an average concentration of 0.01 mg Hg per m# of air in a 
particular industrial operation suspected of causing chronic mercurialism. Mercury in quan- 


tities of 10-8 g was detected, and the method offers possibilities of higher sensitivity and more 


general application. 


INTRODUCTION 


N an investigation of a particular industrial 
operation suspected of causing chronic mer- 

curialism, it was necessary to measure minute 
quantities of mercury vapor under conditions 
where optical measurements'? were not feasible. 
However, it was found that the measurements 
could be made by the use of radioactive mercury* 
as a tracer. Quantities of mercury vapor as small 
as 10-? wg were determined in 2 liters of air. This 
concentration limit was adequate for the problem 
at hand, but the sensitivity could be increased 
by a factor of 10 to 100 if necessary. 

The operation in question consisted of sealing 
small side-arm tubes (see left side of Fig. 1) to 
fluorescent lamp bulbs. 
about 30 mg of mercury. 


The side arm contained 
An operator heated the 
constricted end of the side-arm tube and a cor- 
responding tubulation on the fluorescent bulb 
pressed the softened ends together, and blew 
gently into a mouth piece (not shown) connected 
by a rubber tube to the cotton filter. This 
operation was repeated 50 to 60 times an hour. 
Because of the appreciable 
mercury 


vapor pressure of 


at room temperature, this procedure 


involved the inhalation of some mercury vapor. 
Measurements were made to determine if this 
quantity exceeded the toxic limit of 0.1 mg Hg 
per cubic meter of air.‘ 


'M. Shepherd, S. Schukmann, R. H. Flinn, J. W. 
Hough, and P. A. Neal, J. Research Nat. Bur. Stand. 26, 
357 (1941). 
apa C. Clark and O. Oldenberg, J. Chem. Phys. 9, 786 

). 

3G. Friedlander and C. S. Wu, Phys. Rev. 63, 227 (1943). 

* This toxic limit, established by Shepherd et al. (refer- 
ence 1), is based on a comprehensive experimental study 
and replaces the value of 0.25 mg per cubic meter proposed 
from a review of the literature by C. Goodman, Rev. Sci. 
Inst. 9, 233 (1938). 
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Fic. 1. Diagram of apparatus used in collecting mercury 
from air samples. A known volume of air is drawn through 
the side-arm tube upon connecting the evacuated suction 
flask to the impinger tube. The mercury vapor in the air 
is collected on the refrigerated, amalgamated copper disk. 
This disk is removed, covered with a thin lacquer film, and 
the content of radioactive mercury determined by beta-ray 
counting. The calibration flask was used in an attempt to 
obtain a known concentration of mercury vapor in air. 
This method was not satisfactory. Instead a standard 
solution was used. 


APPARATUS AND METHOD 


Figure 1 is a schematic diagram of the appa- 
ratus designed for these tests.> Radioactive 
mercury was placed in the side arm and air was 
drawn through the tube into an evacuated flask 
of 2.08-1 capacity. Before reaching the flask, the 
mercury-laden air impinged on the amalgamated 
copper disk. This disk, fastened with rubber 
cement to the end of a tube filled with liquid 
oxygen, served to collect the radioactive mercury 
for measurement. 


* A more detailed description of the apparatus and indus- 
trial operations is appearing elsewhere. C. Goodman, J. W. 
Irvine, Jr., and C. F. Horan, J. Ind. Hyg. and Tox., in 
press (1943). 
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After the pressure in the apparatus reached 
atmospheric, the copper disk was removed, 
warmed to room temperature, dried, and covered 
with a thin film (1 mg/cm?) of lacquer. The 
activity of the disk was then measured on a 
thin-window, bell-type counter® with a recording 
counting-rate meter.’ 

An attempt was made to calibrate the appa- 
ratus by collecting the mercury vapor from a 
known volume (112 ml) of air in equilibrium 
with radioactive mercury. This was done by sub- 
stituting a flask containing a small drop of 
mercury for the side-arm tube (Fig. 1). Although 
equilibrium was approached from both sides, the 
observed concentration always bracketed the 
true equilibrium value. These results indicated 
the need for more control to attain equilibrium 
than was considered necessary for this work. For 
purposes of quantifying the activity measure- 
ments in terms of weight of mercury, a standard 
solution of the radioactive mercury was prepared 
by dissolving 0.30 mg in 0.05 M nitric acid. An 
aliquot portion of this solution was evaporated 
on an amalgamated copper disk for measurement. 


PREPARATION OF RADIOACTIVE MERCURY 


Radioactive mercury, Hg!’ (half-life=23 
hours)*, was prepared by bombarding metallic 
gold with deuterons. A yield of 5.6 uc per wamp.h. 
was obtained by probe bombardment with 12- 
Mev deuterons in the M.I.T. cyclotron. 

Five-mil gold foil was silver-soldered to a 
water-cooled probe and bombarded in the 
cyclotron chamber. Before the target was re- 
moved, a piece of 2-mil gold foil was placed over 
the active surface. This foil retained a large 
fraction of the radioactive mercury distilled from 
the target during the heating to melt the solder. 

Both pieces of gold were amalgamated with 85 
mg of inactive mercury carrier and placed in a 
simple vacuum still for purification of the radio- 
active mercury. The still consisted of a 100-mm 
length of 18-mm Pyrex tubing to which a 250-mm 
length of 3-mm capillary tubing was sealed. After 


° A. F. Kip and J. W. Irvine, Jr., Rev. Sci. Inst. (to be 
published, 1943). 


7A. F. Kip, A. G. Bousquet, W. N. Tuttle, and R. D. 
Evans, Rev. Sci. Inst. (to be published, 1943). 
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introduction of the amalgamated gold, the open 
end of the larger tube was sealed and the still 
continuously evacuated with a Hyvac pump. 
Heat was applied to the large tube until the 
mercury distilled into the capillary and the tube 
collapsed. The capillary was then sealed between 
the gold and the condensed mercury. By gentle 
heating the mercury was distilled up the capillary 
in four stages, with a seal-off between each stage. 
It is necessary to carry out these distillations 
without allowing the mercury to boil. In this 
way mechanical transfer of traces of radioactive 
gold, also formed during the bombardment, is 
prevented. The final condensate was collected 
into a small droplet by centrifuging the last sec- 
tion of capillary. 


EXPERIMENTAL RESULTS 


Table I summarizes the experimental results. 
Measurements were made with a 14-mg and a 
70-mg droplet of mercury in the side-arm tube 
and with the 70-mg droplet in the cotton filter, 
samples T-1 to T-9. Room-air blanks, B-1 and 
B-2, and calibration standards, S-1 to S-3, were 
measured. The time recorded in column 3 is the 
number of minutes between the collection of 
successive air samples. The samples were col- 
lected by allowing the evacuated flask to come 
to atmospheric pressure in 35 seconds to 1 
minute. Activity measurements were made in 
duplicate and are corrected for decay of the 
activity and for counter background. The quan- 
tities of mercury are expressed in mg per cubic 
meter (m*) and are calculated from an average 
specific activity of 9690 counts per minute 
(c.p.m.) per wg of mercury. 

In samples T-1 to T-6, the first of each group 
shows a markedly higher activity than the suc- 
ceeding two samples. This decrease is probably 
caused by the formation of an oxide film on the 
mercury while standing undisturbed and the 
lack of time for mercury vapor to diffuse into the 
straight portion of the side-arm tube after this 
has been freed of the vapor that accumulated 
during transfer of the droplets. With the droplet 
in the direct air stream, T-7 to T-9, the successive 
samples are more consistent. 

The high blank values of the room air probably 
resulted from the fact that the apparatus was 
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TABLE I. Summary of numerical data from mercury tests. 


Counts per 
min. from 
min.* 2.08 

of air 


Time 
Mercury Sample in Hg in mg/m? 
data No. i 


Observed Average 


14 mg in 655 
side arm 715 0.037 
of tube T 65 0.0034 

61 0.0031 
14 | 0.0007 
26 =|: 0.0013 


0.034 0.036 


0.0033 


0.0010 


70 mg in 7 ; 367 
side arm 358 | 0.018 
of tube T-§ 7 185 0.0095 

186 | 0.0095 
80 0.0041 
105 | 0.0054 


0.019 0.019 


0.0095 


0.0048 


70 mg in T . 123 
cotton 318 | 0.016 
filter T-8 280 0.014 

284 | 0.014 
204 | 0.011 
212 | 0.011 


0.0063 0.011 


0.014 
0.011 


Blank 32 =| 0.0016 
36 0.0019 
Blank - 99 0.0051 


101 0.0052 


0.0018 


} 


0.0052 


| Counts per min. 


per 10°* g Hg 

| Observed 
4600 
5330 
5850 


Average 


9,200 
10,660 
9,160 


5x10~* 
mg Hg 

6.410 
mg Hg 

13x 10~* — 
mg Hg 


12,440 9,750 


* This is the length of time the metallic mercury was in contact with 
the air in which the mercury vapor was measured. 


set up in the radiochemical laboratory where the 
mercury was purified. Although the separation 


was made under a good fume hood, the mercury 
was divided and weighed in the room and suf- 
ficient mercury vaporized to contaminate the air. 

Measurements of the activity on the outside 
of the liquid-oxygen tube after the collection of 
a very active sample (10,000 c.p.m.) from warm 
mercury were made. These results indicated that 
the removal of mercury vapor from the air by 
this apparatus was more than 99 percent com- 
plete. 

The highest concentration of mercury vapor 
obtained from the side-arm tube was only 36 
percent of the toxic limit of 0.1 mg per m* set 
by Shepherd et al.1 The average was about 10 
percent of this value. From these results it may 
be concluded with some certainty that under the 
operating conditions in question toxic quantities 
of mercury vapor could not be inhaled. 


OTHER APPLICATIONS 


The experiments reported in this paper rep- 
resent a highly specialized application of the use 
of radioactive mercury. However, the extreme 
sensitivity (<10-*§ g Hg) of the method suggests 
other applications where minute quantities of 
mercury are of interest. For instance, absorption 
of mercury through the skin, lungs, stomach, or 
teeth could be accurately measured. The rate of 
elimination of mercury from the body and 
factors which control that rate are subject to 
investigation by the radioactive tracer method. 
Other applications undoubtedly suggest them- 
selves to the reader. 
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